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An accurate quantum three-body calculation is performed for the new type of big-bang 
nucleosynthesis (BBN) reactions that are catalyzed by a hypothetical long-lived negatively 
charged, massive leptonic particle (called X~) such as the supersymmetric (SUSY) particle 
stau, the scalar partner of the tau lepton. It is known that if the X~ particle has a lifetime of 
Tx ^ 10'^ s, it can capture a light element previously synthesized in standard BBN and form 
a Coulombic bound state, for example, (^BeX~) at temperature Tg <^ 0.4 (in units of 10^ K), 
{aX~) at Tg <, 0.1 and (pX~) at Tg <^ 0.01. The bound state, an exotic atom, is expected 
to induce the following reactions in which X~ acts as a catalyst: i) a-transfer reactions such 
as [aX~)-\-d — > ^\A-\-X, ii) radiative capture reactions such as (^BeX~)+p (*BX~)-|-7, 
iii) three-body breakup reactions such as C\AX~) + p ^ a + a-\- X~ , iv) charge-exchange 
reactions such as {pX') -I- q — > (aX^) + p and v) neutron induced reactions such as 
{^BeX-) + n ^ ^Be-I-X". In recent papers it has been claimed that some of these X - 
catalyzed reactions have significantly large cross sections so that the inclusion of the reactions 
into the BBN network calculation can markedly change the abundances of some elements, 
giving not only a solution to the problem (the calculated underproduction of ®Li by 

a factor of ~ 1000 and overproduction of ^Li-|-^Be by a factor of ~ 3) but also a constraint 
on the lifetime and primordial abundance of the elementary particle X~ . However, most of 
these calculations of the reaction cross sections in the literature were performed assuming 
too naive models or approximations that are unsuitable for these complicated low-energy 
nuclear reactions. We use a high-accuracy few-body calculation method developed by the 
authors and provide precise cross sections and rates of these catalyzed BBN reactions for 
use in the BBN network calculation. 

§1. Introduction 

The accurate theoretical prediction of nuclear reaction rates is of essential impor- 
tance in the study of big-bang nucleosynthesis (BBN) when particular reaction rates 
are very difficult or even impossible to estimate experimentally. The most typical ex- 
ample of such reactions is the production and destruction of light elements catalyzed 
by a hypothetical long-lived negatively charged massive 100 GeV) leptonic parti- 
cle, here denoted as X~ ; a candidate for this particle is the supersymmetric (SUSY) 
counterpart of the tau lepton (r), i.e., the stau (f).*) Recently, BBN involving these 
X^-catalyzed nuclear reactions has been extensively studied from the viewpoint of 
particle-physics catalysis in BBN,"*^-*"^^^ which may be briefly reviewed as follows (and 
precisely reviewed in Refs. 12) and 14) and references therein): 

i) If the X~ particle has a lifetime on the order of 10^ s or longer, it can capture 
a light nucleus, A, previously synthesized during BBN, and produce a bound state 



*' The stau is one of the particles expected to be discovered at the CERN Large Hadron Col- 
lider.is) 
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denoted as {AX~), which is a type of exotic atom with the nucleus ^ in a Coulombic 
orbit around the massive X~ at the center.*^ 

ii) When the cosmic temperature cools to Tg ~ 0.4 (in units of 10^ K) after the 
standard BBN (Tg S> 0.5), X~ captures '^Be to form (''BeX^), which is subsequently 
destroyed by the (''BeX") + p —>■ (®BX~) + 7 reaction^) followed by the /3-decay of 
{^BX-) to a + a + X-. 

iii) At Tg ~ 0.3, X~ captures ^Li to generate (^LiX~), which is subsequently 
destroyed by the (^LiX~) + p ^ a + a + X~ reaction. 

iv) At Tg ~ 0.1, X^ particles capture a particles to form {aX~) in abundance 
which triggers the reaction {aX~) + d — ^ ^Li + X" to produce an enormous amount 
of ^Li, as originally proposed by Pospelov^^ and confirmed by Hamaguchi et al.^^ 
including two of the present authors (M. K. and Y. K.). 

v) At Tg ~ 0.01, although X^ particles start capturing protons to form {pX~), its 
abundance remains low due to the charge-exchange reaction {pX~)+a {aX~)+p, 
which occurs before {pX~) begins to interact with other nuclei to form heavier 
elements. 

The contribution of these reactions is expected to modify standard BBN (SBBN 
in short) and to solve the calculated underproduction (by a factor of ~1000) of 

the primordial abundance of ^Li and the overproduction (by a factor of ~3) of 
^Li,^^^ and, at the same time, to provide information for constraining the lifetime 
and primordial abundance of the elementary particle X~ . However, full quantum 
mechanical studies of the reactions have not been performed on the basis of nuclear 
reaction theory except for the work of Rcf. 5) on the production of ^Li. 

The purpose of the present paper is to perform a precise quantum three-body 
calculation of the cross sections of many types of catalyzed BBN (CBBN in short) 
reactions including those mentioned above in i) - v) and to provide their reaction 
rates for use in the BBN network calculation. The adopted theoretical method 
and numerical technique have been developed by the present authors, reviewed in 
Ref. 17), and are well established in the field of few-body atomic and nuclear systems. 

The CBBN reactions studied in this paper are classified into four types: 

1) Transfer reactions: 

{AX-)+a^ B + X-, (1-1) 

where A is picked up by an incoming nucleus a to produce a nucleus B{= A + a), 

2) Radiative capture reactions: 

{AX-) + a^ {BX-) + -f, (1-2) 

3) Three-body breakup reactions: 

{AX-) + a^bi + b2 + X-, (1-3) 

4) Charge-exchange reactions: 

{AX-) + a^ {aX-) + A. (1-4) 

*^ For instance, the binding energies of (^BeX~) and {aX~) are ~1.3 and ~0.34 MeV, and the 
rms radii are ~3.6 and ~7 ftn, respectively. 
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We solve the Schrodinger equations of CBBN reactions (1-1)— (1-4), exphcitly taking 
the three-body degree of freedom into account, and derive their cross sections and 
reaction rates at incoming energies below 200 keV. 

Throughout this paper, it is not necessary to identify X~ with any particular 
particle such as the SUSY particle stau, although this is assumed, for instance, in 
Refs. 4), 5) and 14). The property of X~ we assume here is that it has a charge of 
— e, a mass mx ^ 100 GeV and a lifetime tx ^ 10^ s and that it interacts with light 
elements via the Coulomb force only. 

The structure of the present paper is as follows. We investigate X~ -catalyzed 
a-transfer reactions that produce ^Li, ^Li and ^Be in §2, resonant and nonresonant 
radiative capture reactions that destroy ^Be in §3, three-body breakup reactions that 
destroy ^Li and ^Li in §4, the charge-exchange reactions in §5 and the possibility of 
the X~-catalyzed primordial production of ^Be in §6. A summary is given in §7. 



§2. X -catalyzed a-transfer reactions 



2.1. Necessity of three-body calculation 
2.1.1. Production of ^Li 

One of the puzzles in SBBN is that the BBN prediction of the ^Li abundance 

is a factor of ~1000 smaller than the observed data.^^^ The too small prediction is 
mostly because the radiative capture reaction that produces ^Li in SBBN, 

a-\-d^ ^Li + 7 (> 1.47 MeV), (2-1) 

is heavily El-hindered with a very small astrophysical S-factor, S'-y ~ 2 x 10""^ 
keV b.^^) To obtain a much greater production of ^Li, Pospelov considered,^) for the 
first time, a new type of CBBN reaction, 

(aX") ^Li + X" + 1.13 MeV, (2-2) 

in which the bound state*) {aX~) promotes the synthesis of a and d to form ^Li 
with no photon emission. He claimed that the S'-factor of the -catalyzed reaction 
(2-2), say Sx, is enhanced by 8 orders of magnitude {Sx ~ 3 x 10^ keV b) compared 
with Sj and that the abundance of ^Li is greatly increased so that the ^Li problem 
can be solved under a suitable restriction on the abundance and lifetime of X~ . 

His estimation of Sx, however, was based on a too naive comparison between 
(2-1) and (2-2). Namely, in an analogy to the real photon in (2-1) with the wavelength 
A-y ~ 130 fm, Pospelov introduced a virtual photon in (2-2) with a characteristic 
wavelength Ax on the order of the Bohr radius (= 3.63 fm) of the atom (aX~) and 
assumed the scaling relation**) 



*^ Hereafter, when {AX~) is in the Is ground state, the suffix Is of {AX~)is is omitted. 
**^ A more detailed expression is given by Eq. (4) in Ref. 1) but the scaling relation is essentially 
described by (2-3). 
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where the photon multipolarity i is assumed to be £ = 2 (namely, E2), not only for 
the real photon*^ but also for the virtual one. In this case, the enhancement ratio 
Sx/Sy amounts to ~ 10^. 

However, the scaling relation (2-3), which compares the wavelengths of the real 
and virtual photons, is not suitable for discussing reaction (2-2) because the relation 
does not take into account the nuclear interaction, which is the most important 
factor in the process of transferring the a particle from to the deuteron to 

form ^Li. Therefore, two of the authors (M. K. and Y. K.) and their colleagues^-* 
performed a fully quantum three-body calculation of the a + d + X~ system and 
calculated the /S-factor of reaction (2-2), Sx{E), as a function of the incident energy 
E at E = 10 — 120 keV. The calculated S'-factor at the Gamow peak energy is 38 
keV b (Sx/Sj ~ 10^, nearly one order of magnitude smaller than that in Rcf. 1)). 
Note that Sx of the CBBN reaction (2-2) is of a similar order of magnitude to the S- 
f actors of nonresonant photonless SBBN reactions caused by the nuclear interaction; 
for example, compare Sx with the observed values of ~ 60 keV b for d+d t+p and 
'--^ 5 X 10~^ keV b for the El-forbidden process d+d ^ a + 'j. The large enhancement 
of the CBBN/SBBN ratio (~ lO'^) is simply because the SBBN reaction (2-1) is 
heavily El-hindered. 

2.1.2. Production of ^Li and ^Be 

In the BBN network calculation with the CMB-based baryon-to-photon ratio 
(?7b = 6.0 X 10^^*^)."'^^-' the abundance of '^Bc is approximately one order of magnitude 
larger than that of ' Li. Since ^Be eventually decays to ''Li by electron capture, the 
sum of the abundances of both ^Li and ^Be is regarded as the abundance of ^Li in 
the ^Li-overproduction problem. To examine the effect of CBBN on the abundance 
of ^Li-h^Be, Cyburt et al.*^ studied the a-transfer reactions 

{aX-)+t^ ^Li + X- +2.13MeV, (2-4) 
(aX-)+^He^ ^Be + X- + 1.25MeV. (2-5) 

which produce ^Li and ^Be at Tg ~ 0.1. Using the scaling relation proposed in 
Ref. 1), they claimed that the S'-factors of the above CBBN reactions are 5 orders 
of magnitude greater than those of the El radiative capture reactions 

a + t^ ^Li + 7(> 2.47MeV) , (2-6) 
a + ^He^ ^Be + 7(> 1.59MeV), (2-7) 

which are the most effective source for producing ^Li and ^Be in SBBN. 

However, it is clear that the enhancement ratio CBBN/SBBN in S'-factors will 
be only on the order of 10^ — 10^ contrary to the case of ^Li production (~ 10^). The 
reason for this is as follows: the SBBN S'-factors of (2-6) and (2-7) at astrophysical 

*^ This assumption oi i = 2 for the real photon in (2-1) may not be justified at astrophysical 
energies (no experimental data there). According to the so far most precise six-body calculation^"' 
of the El-hindered radiative capture reaction (2-1), the S-factor of the E2 transition decreases much 
more rapidly than that of the El transition as the energy decreases, and eventually the E1/E2 ratio 
becomes ~ 10 at £ ~ (see Figs. 7 and 8 in Ref. 20)). 
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energies are respectively ~ 0.1 and ~ 0.5 keV b/^-* which are 5 orders of magnitude 
larger than that of the El-hindered reaction (2-1). Moreover, the CBBN S'-factors of 
(2-4) and (2-5) should be significantly reduced compared with that of (2-2), because 
a rearrangement of the angular momenta among the three particles is needed when 
going from the entrance channel to the exit channel; note that a significant difference 
between (2-2) and (2-4)— (2-5) is that the relative motion between d and a in ^Li is 
that of s-wave whereas that between t(^He) and a in ^Li (^Be) is that of p-wave. 
This consideration of the CBBN/SBBN ratio will be confirmed by the quantum 
three-body calculation of reactions (2-4) and (2-5) in this section. 

2.2. Three-body Schrodinger equation and reaction rate 

Following Refs. 5) and 17), we briefly explain the three-body calculation method 
used to investigate the above X~-catalyzed transfer reactions. We emphasize that 
the method described here is an exact method for calculating the cross section of 
the low-energy transfer reactions within the three-body model of the total system 
concerned. The present authors have previously performed the same types of cal- 
culations in the study of muon transfer reactions^^)'^-*^)"^^) in muon-catalyzed fusion 
cycles. 

As shown in Fig. 1, we consider all three sets of the Jacobi coordinates (vc, Rc), 
c = 1, 2 and 3, to completely treat the three-body degrees of freedom of the system 
A-\- a -\- X~. Here, we take A = a, and a = d,t and ^He for reactions (2-2), (2-4) 
and (2-5), respectively. 



A A A 




Fig. 1. Three sets of Jacobi coordinates in the A+a+X~ system. The entrance channel {AX~)+a is 
described using the coordinate system c = 1, and the transfer channels (Aa)+X~ and {aX~)+A 
are described using c = 2 and c = 3, respect vely. 

The Schrodinger equation for the total wave function with total energy £^tot, 
angular momentum J and ^-component M is given by 

{H - Etot)^jM = (2-8) 

with the Hamiltonian 

^ = -^^-c - ^^R. + VA-x{r^) + VA-a{r2) + Va-xir^). (2-9) 
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As far as we use the reduced masses {nic and Mc) associated with the coordinates 
{Vc and Rc), every choice of c in the kinetic term is equivalent. The notation of 
the potentials is self-evident; an explicit form will be given below. The Schrodinger 
equation is solved under the scattering boundary condition imposed appropriately 
on the relevant reaction. 

For all the CBBN reactions studied in this paper, only the resonant radiative 
capture reaction (1-2) requires serious consideration of the intrinsic spins of the 
particles as will be discussed in §3. In this section, for the nonresonant a-transfer 
reactions, the spin of the incoming particle a (= d, i,^He) is neglected. Therefore, in 
the exit channel of (2-4) and (2-5), the 3/2~(p-wave) ground state and the l/2~{p- 
wave) first excited state of ^Li and ^Be are degenerated energetically to their weighted 
mean. The validity of this approximation will be discussed in §2.6. 

Firstly, we construct the wave function of the Is ground state of (AX^ ) in the 
entrance channel on the coordinate system c = 1 and that of the ground state of the 
nucleus B on c = 2. We denote the wave function by <?^[^^^(^c) and the eigenenergy 

by elf; explicitly, li = for (aX'), h = for ^Li and h = 1 for '^Li (^Be). These 
quantities are obtained by solving 



[-^V^^ + K(rc)-e!f]0ll(^c) = O, (c=l,2) (2-10) 



where the potential Vc denotes Va-x and Va-q for c = 1 and 2, respectively. 

Using XlImS^c^ ^ (-^c)^LcMc(-Kc), we denote the scattering wave function 
along the coordinate Rc with angular momentum Lc and z-component Mc (c = 1,2). 
The center-of-mass (cm) scattering energy associated with the coordinate Rc, say 
Ec, is introduced as 

Ec = Etot - e^f = f?kl/2Mc, (c = 1, 2) (2-11) 

together with the corresponding wave number kc- 

The total three-body wave function iP'jm, which describes the transfer reaction 
(2-2) in the most sophisticated manner, is written as"'^^-'''^^-' 

^JM = 0S (ri) x%{R,) + [^ff (r,) xg (^2)] + • (2-12) 

Here, the first and second terms represent the two open channels, {AX~) -f a on 
c = 1 and B + X~ on c = 2, respectively; there are no other open channels in the 
energy range (Ei < 150 keV) that we are interested in. In the first term, the radial 

part Xj^(i?i) has incoming and outgoing amplitudes, and Xl2(-^2) in the second 
term has the outgoing amplitude only. These should satisfy the boundary condition 

hm RcX^^liRc) = u^£-Jikc,Rc)Sci-.I^S^^^u^^Jikc,Rc), (c=l,2) (2-13) 

where Li = J, and u^^\kcRc){= G L{kcRc)^'iFL{kcRc)) are the asymptotic outgoing 
and incoming Coulomb wave functions. S^^^ is the S'-matrix for the transition from 
the incoming channel (c = 1) to the outgoing channel c and Vc is the velocity. 
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The third term in (2-12), represents all the closed (virtually excited) 

channels in the energy range in this work; in other words, this term is responsible for 
all the asymptotically vanishing amplitudes due to the three-body degrees of freedom 
that are not included in the first two scattering terms.*) Since ^j^m^'^^^ vanishes 
asymptotically, it is reasonable and useful to expand it in terms of a complete set of 
L^-integrable three-body basis functions, {^jm,u'iI^ = 1 — t'max}, that are spanned 
in a finite spatial region (see §2.4): 



max 



^(Closed) ^Vftj.^.M,.. (2-14) 



Equations for Xj)(i2i), Xi^i^"^) coefficients bji, are given by the i^max + 2 

simultaneous equations, ^^^'^^^ 

( cP^^!,{ri)YjM{Ri) I H - Etot I ^JM = 0, (2-15) 

( (^2) ® Yl, (R^)] ^JH- E,,, I ),^, = 0' (2- 16) 

and 

( <PjM, v\H-E^^t\^JM)=^- (v = \- v^^) (2-17) 

Here, ( ) ^ ^ denotes integration over and Rc- 

Since ^jm, v are constructed so as to diagonalize the three-body Hamiltonian as 
(cf. §2.4) 

{^jM,v\ H \^jMy) = Ej^6,,„>, (z/, u' = 1- fmax) (2-18) 
the coefficients bj,y can be written, from Eq. (2-17), as 



bj, = ^(^jM,. I H - Etot I 4oVi)x?i(i2i) + ® xS(i«2) 

(Z^=l-I^ma.) (2-19) 

Using Eqs. (2-15), (2-16) and (2-19), we reach two coupled integro-differential equa- 
tions for x^j\Ri) and x[j\R2), which are not written here (see §8 of Ref. 17) for 
the equations). 

Finally, the integro-differential equations are solved using both the direct numeri- 
cal method (the finite-difference method) and the Kohn-type variational method. ^^^'^^^ 
We obtained the same result; this demonstrates the high accuracy of our three-body 
calculation. 

Using the S'-matrix elements obtained above, the cross section of the rearrange- 
ment process is expressed as 

00 

'^(^) = pE(2^ + 1)1^1-21'- (2-20) 

'^i j=o 



*^ This method for low-energy three-body reactions has already been used^^''^^'"^^^ in the study 
of the muon transfer reaction (d/i)is -h t — > (t/i)is -h d -|- 48 eV. The ^j'^m^"'^^ term was found to play 
a very important role; if the term is omitted (namely, if the two-channel coupled calculation is 
performed with c = 1, 2), the calculated low-energy (0.001 — 100 eV) cross section of the reaction 
becomes ~ 30 times larger than that obtained by the full three-body calculation. 
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where we introduce a simplified notation, E = E\, for the energy of the entrance 
channel. The astrophysical S'-factor is derived from 



a{E) = S{E)exp{-27rr]{E))/E. 



(2-21) 



Here, cxp{—2TTrj{E)) is the Coulomb barrier penetration probability, where r]{E) is 
the Sommerfeld parameter of the entrance channel {AX~) + a, defined as 



ri{E) = ZaZAxe^/hvi, 



(2-22) 



where ZqC and Zaxg are the charges of a and {AX ), respectively. 

The reaction rate {(tv) at temperature T is expressed as (cf. Eq. (4-44) of 
Ref. 26)) 



{av) 



where k is the Boltzmann constant. If S{E) is simulated by a linear function of E 
around the Gamow peak energy Eq as 



S{E)c^S{Eo)+l^^^^ (E-Eo), 



(2-24) 



Eo = 122.0 {Z^Zlxl^)^ keV, (2-25) 
the reaction rate Na {o'v) is expressed, using Eqs. (4-56) and (4-75) of Ref. 26), as 



Na{<tv) = 7.82 X 10' 



X Tg ^ exp 



6 / ZoZax ^ ^ 



S{Eo) + 71.8 ( ^ j Tg 



- 4.248 (z2zl;,/x)iT9 



cm^s'^mor^ (2-26) 



where S{Eo) and {dS /dE)Eo are given in units of keVb and b, respectively, Na is 
Avogadro's number and is Mi in units of amu. As will be seen for any nonresonant 
reaction in this paper, we can regard {dS/dE)Eo ~ constant = a: 



S{E) = S{Eo) + a{E-Eo) = S{0) + aE, 



(2-27) 



where S{0) and a are given in units of keV b and b, respectively, and E is in keV. 
We then have 



Na{ctv) 



.82x10^ (M^y siO)T;hxp 



4.248(Z2Z1^/X)ir9 3 



1 + 



a 



5(0) 



122.0 (Z„^Zi^/x) 3 Tg^ +71.8 Tg 



cm^s~Vor^. (2-28) 
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2.3. Nuclear and Coulomb potentials 

It is essential in the three-body calculation to employ appropriate interactions 
among the three particles. Here, we determine the potentials VA-x{ri), VA-a{r2) and 
Va-x{f3) in (2-9), where A = a and a = d,t and '''He. 

Wc assume the Gaussian charge distribution Ze{TTb'^)^'^^'^e^^^^'^^ for a, d, t and 
^He; here, Ze is the charge and b = -^2/3 tq, where tq is the observed rms charge 
radius, which is given by ro = 1.67, 2.14, 1.70 and 1.95 fm, respectively. 

The Coulomb potential between A and X~ is given by 

VA-xir) = -ZAe^'-^^, (2.29) 
r 

where erf(x) = e'^'dt is the error function; and the Coulomb potential is 

similarly defined for Va-x{r)- The energy of the {aX~)is state is s'gs = —337.3 keV 
(—347.6 kcV) and the rms radius is 6.84 fm (6.69 fm) for mx = 100 GeV {mx oo). 

The potential VA-a{^) is a sum of the nuclear potential, V^^(r), and the Coulomb 
potential, V^_g^{r). The latter is given by 



V^-a{r) = ZAZa e^ ^— . (2-30) 

The nuclear potential V^_^{r) is assumed to have a two-range Gaussian shape as 

ylair) = VI e-^''/"!)' + V2 e-('"/»2)'. (2-31) 

i) a-d potential 

Wc take ai = 0.9 fm, vi = 500.0 MeV, 02 = 2.0 fm and V2 = -64.06 MeV.^) The 
first term, a repulsive core, is introduced to simulate the Pauli exclusion principle 
that nucleons in the incoming deuteron should not occupy the nucleon s-orbit in 
the a particle during the reaction process (for this role of the Pauli principle,*-* see, 
for example, Ref. 27)). The parameters are determined so that the solution to the 
Schrodinger equation (2-10) for c = 2 reproduces observed values of both the energy 

Sgs^ = —1.474 McV and the rms charge radius 2.54 fm^^^ of the ground state of ^Li. 
Furthermore, the charge density of ^Li reproduces the observed charge form factor 
of the electron scattering from ^Li (see Fig. 2 of Ref. 5)). Simultaneously, the use 
of the potential Va-d{r2) explains the low-energy s-wave phase shifts of the a + d 
scattering (see Fig. 2 of Ref. 5)). 

ii) a-t potential 

The nuclear potential between a and t is assumed to be parity dependent; ai = 1.0 
fm, vi = 500.0 MeV, 02 = 2.7 fm and V2 = —46.22 MeV for odd angular-momentum 



*^ Use of the sophisticated orthogonality-condition model (OCM)^^' for three-body systems is 
not necessary in the present case, since the Pauli principle only applies between A and a among the 
three particles, and we are not treating their compact bound states. However, the introduction of the 
inner repulsive potential is useful in this type of three-body scattering calculation to automatically 
prevent the unphysical Pauli forbidden amplitude from entering the total wave function. 
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states and ai = 1.5 fm, vi = 500.0 MeV, 02 = 2.4 fm and V2 = -7.0 MeV for 
even angular-momentum states. The repulsive term is introduced so as to prevent 

the Pauli forbidden states from entering the total wave function. The use of the 

(2) 

potential reproduces the observed energy Sgs = —2.308 MeV (weighted average for 
the 3/2~ ground state and the l/2~ excited state) and the rms charge radius 2.43 
fm^^) of the ground state of ^Li as well as the observed values of the low-energy a-t 
scattering phase shifts for the p-wave (average of the and pi/2 states) and the 
s-wave. 

iii) a-^He potential 

The nuclear a-^He potential is assumed to have the same shape and parameters as 
the a-t potential except that V2 = —44.84 MeV for odd states, which reproduces the 
observed weighted mean energy {e^s = —1.444 MeV) of the ground state and the 
l/2~ excited state. 

2.4. Three-body basis functions 

The L^-integrable three-body basis functions {<P,jMy, ly = 1 — i^max} used in 
(2-14) to expand li^Sr"^^ are introduced as follows: '^jm v are written as a sum 
of the component functions in the Jacobi coordinate sets c = 1 — 3 (Fig. 1), 

#JM,. = #?i,,(ri,iii) +^S^^,(r2,i^2) +^il,,(r3,i23) - (2-32) 
Each component is expanded in terms of the Gaussian basis functions of and Re- 

^Sl.(^c,i^c)= E 4tne^c,iV.Le [C/c(^c)V'^.L.W]^M'(c=l-3) (2-33) 

where the Gaussian ranges are postulated to lie in a geometric progression: 

<^nlrrXr) = r' e'^^l^^^^ YUf), Tn = n a^-\ (n = 1 - n^ax) (2-34) 
i^NLMiR) = R"" e-WRN? Ylm{R), Rn = Ri A^'-K {N = 1 - N^^) (2-35) 

The basis functions chosen in this way arc suitable for describing both short-range 
correlations (mainly due to nuclear interactions) and long-range asymptotic behavior 
simultaneously, and therefore they are efficient for describing any three-body configu- 
ration (the closed-channel contribution) in the interaction region of the intermediate 
stage of reactions. The coefficients A^fl ^ ^ ^ in (2-33) and the eigenenergies E'j 
of ^jM,u are determined by diagonalizing the three-body Hamiltonian H as (2-18). 

This method in which the total Hamiltonian is diagonalized using the precise 
Gaussian basis functions has successfully been applied to the study of various types 
of three- and four-body systems in physics,*^ which is reviewed in Ref. 17). 

i) Basis set for a -\- d-\- X~ 

The most precise three-body calculation among the applications is the determination^'^^'^"^ 

of the antiproton mass''^-' by analyzing CERN's laser spectroscopic data on the antiprotonic helium 
atom, He++ + The energies of the three-body atom with J ~ 35 were calculated with an 

accuracy of 10 significant figures. 
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In the calculation for J = 0, we took lc = Lc = 0, 1, 2 and nmax = ^max = 15 for c 
1 — 3. The total number of the three-body Gaussian basis [(j)^ ^ {tc) i {R, 



JM 

used to construct the cigcnfunctions {(PjM,u} amounts to J^max = 2025, which was 
found to be sufficiently large for the present calculation. For the Gaussian ranges, we 
took {ri,rn^^, Ri, RNma.^} = {0.5, 15.0, 1.0, 40.0 fm}, which are sufficiently precise 
for the present purpose. The expansion (2-14) converges quickly with increasing u, 
and 1/ < 100 {Eji, < 1 MeV above the {aX~) — d threshold) is sufficient for the 
conversion. Basis sets for J > 1 are not shown since the contribution to the cross 
section from J > 1 is 3 orders of magnitude smaller than that from J = 0. 

ii) Basis set for a + t (^He) + X~ 
Since Z2 = 1 in the exit channel of (2-4) and (2-5), the basis [^^^^^('"c) V'^^l^(-Rc)] 
for J = 1 with {(^1 = 0, -Li = 1), {I2 = 1. ^2 =0)} is the most important (sec Figs. 3 
and 4). The next most important set for J = 1 is {{li = l,Li = 0), {I2 = 1,1^2 = 
0)}. We take Gaussian ranges {ri,rn^^, Ri, RN^a^} = {0.4, 15.0, 0.8, 30.0 fm} with 
nmax = -/Vmax = 12. The Contribution of J = to the S-f actor is several times 
smaller than that of J = 1; we take the set {Ic = Lc = 0, 1) for c = 1 — 3. 

2.5. Result for ^Li production 

In this subsection we study the ^Li production reaction (2-2) that is the most 
important CBBN reaction. A three-body calculation of (2-2) was reported in Ref. 5) 
by Hamaguchi et ai, and therefore we simply recapitulate it here. The calculated 
astrophysical S'-factor S{E) is shown in Fig. 2 together with the Gamow peak at 
£^0 = 33 keV for Tg = 0.1 {kT = 8.62 keV) around which {aX') is formed. The 
contribution from the s-wave incoming channel with (li = Li = 0) ^ (^2 = L2 = 0) 
is dominant and that from the p-wave incoming channel with (/i = 0, Li = 1) 
(Z2 = 0, L2 = 1) is 3 orders of magnitude smaller. A significant effect of the closed- 
channel amplitude iP'j'j'"^'^*^'' in (2-12), which represents the contribution of the three- 
body degree of freedom in the interaction region, can be seen as follows. If the 
term ^jf^^^^^ is omitted in the three-body calculation, S{E) becomes nearly 3 times 
smaller than that in Fig. 2. 

The CBBN 5-factor in Fig. 2 is enhanced by a factor of ~ 10^ compared with that 
of the SBBN reaction (2-1). This confirms the large enhancement (~ 10^) pointed 
out by Pospelov,^-* although his scaling relation for the estimation is too naive. Note 
that the S'-factor in Fig. 2 is of a similar magnitude to that in typical nonresonant 
photonless SBBN reactions and that the large enhancement ratio originates simply 
from the fact that the SBBN reaction (2-1) is heavily El-hindered. 

The energy dependence of S{E) in Fig. 2 is approximated by (2-27) with S'(O) = 
44.6 keV b and a = —0.18 b. The reaction rate is then given, using (2-28), as*) 

2 12 

NA{crv) = 2.78 x 10^ exp (-5.33T9"3 )(1 - 0.62r93 - 0.29r9) cm^ s'^ mol"^ 

(2-36) 



*' The expression (2-36) is slightly different from (4-3) in Ref. 5). The latter is for the vicinity 
of kT = 10 keV (Tg = 0.116), while the former is vahd for a wider range of Tg ^ 0.2. 
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for Tg < 0.2. On the other hand, since the SBBN reaction rate of (2-1) is given^^^ as 

NA{(7v)sBBf^ = 14.8r9"^exp(-7.44r9"^)(l+6.57r9 + ---) cm=^s-Vor\ (2-37) 

the CBBN/SBBN ratio of the reaction rates is ~ 10^ at Tg = 0.1. This increase of 
the ratio from that of the S-factors (~ 10^) is due to the reduction of the CBBN 
Coulomb barrier experienced by the incoming particle (deuteron) in (2-21). 

60 



1 40 
CO 

20 



I I I I I 
(a X") + d - \\ + X" 




50 

E [keV] 



100 



Fig. 2. The astrophysical 5-factor of the CBBN reaction (2-2) obtained by the three-body calcula- 
tion (solid line). The dotted curve (in arbitrary units) illustrates the Gamow peak for Tg = 0.1 
(fcT = 8.62 keV) with the maximum at Eo = 33 keV. This figure is taken from Ref. 5). 



It is desirable to quantitatively examine how sensitive the calculated cross sec- 
tion (S'-factor) is to the choice of the a-d potential parameter set that reproduces 
the empirical a-d binding energy and the measured rms charge radius of ^Li. Wc 
examined this sensitivity by adopting another parameter set {ri = 0.9 fm, Vi = 400 
MeV, r2 = 2.5 fm, V2 = —35.04 MeV}, which is markedly different from the set 
shown in §2.4. However, we found that the new result for the S'-factor differs by only 
5 — 6% from that in Fig. 2. One can expect a similar result for other nonresonant 
reactions studied in the present paper, although such a test is not repeated there. 

So far, the mass of the X~ particle, mx, has been assumed to be 100 GeV. 
However, we found that the calculated S'-factor differs by only ~ 0.5 % between 
mx = 100 GeV and mx 00. Therefore, we take mx = 100 GeV throughout the 
present paper except for §3, where a resonant reaction is studied. 

Using this large CBBN rate (2-36), Hamaguchi et al.^^ solved the evolution 
equation for the ^Li abundance after SBBN has frozen out. This calculation was 
performed for various sets of values of the assumed lifetime tx and the initial number 
density Ux- of the X~ particle together with the calculation of the number density 
of the bound state {aX~) as a function of the temperature (cf. Fig. 5 of Ref. 5)). 
For the limiting case of a long lifetime {tx 3> 1000 s), they obtained ne^/nB — 
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3.7 X 10~^ nx-/nB, where neLj and hb are the number densities of ^Li and baryons, 
respectively. Therefore, the observational upper bound for the ^Li abundance ^Li < 
6.1 X 10~^^ (2(t)^^^ leads to a remarkable bound for the abundance, nx- /n-Q < 
1.6 X 10~^. The impHcation of this result for particle physics is discussed in Ref. 5). 

2.6. Result for ^lA and ^Be production 

The calculated astrophysical S'-factor of the ^Li production reaction (2-4) is 
shown in Fig. 3. The magnitude of S{Eq) at the Gamow peak energy is approxi- 
mately one order smaller than that of the ^Li production reaction in Fig. 2. This is 
because the ground state of ^Li has s-wave (Z2 = 0) angular momentum of the d-a 
configuration whereas that of '''Li has p-wave {I2 = 1) angular momentum between 
t and a. The incident wave with = Li = must cause the angular-momentum 
rearrangement to Z2 = ^2 = 1 in the exit channel, and therefore the most effective 
partial wave is that with Zi = 0, Li = 1 in the incident channel and Z2 = 1) -^2 = in 



> 

^ A 



(a X") + 1 



'Li + X" 



p-wave 



s-wave 



50 
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Fig. 3. The S-factor of the CBBN reaction (2-4) obtained by the three-body calculation (the solid 
line). The s- and p-wave contributions are shown individually. The dotted curve (in arbitrary 
units) illustrates the Gamow peak for T9 = 0.1 {kT = 8.6 keV) with Eo = 38 keV. 



the exit channel. This can be seen in Fig. 3; the contribution of the p-wave (Li = 1) 
in the incident channel is much larger than that of the s-wave (Li = 0). 

The large effect of the closed-channel amplitude ^f^^^^^ in (2-12), which repre- 
sents the contribution of the three-body degree of freedom in the interaction region 
can be seen as follows. If the term ^j^^^^'^^ is omitted in the three-body calcula- 
tion, S{E) becomes nearly 50 times smaller than that in Fig. 3. The effect changes 
markedly from the case of ^Li production in §2.5 (namely, from 3 times to 50 times). 
This is due to the fact that the angular-momentum transfer between the entrance 
and exit channels is difficult at low energies much lower than the Coulomb barrier, 
and therefore the transfer is strongly mediated by the degree of three-body distortion 
in the internal region where the reaction takes place. 
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In the above calculation, the spin of t is neglected and therefore the 3/2~ ground 
state and l/2~ excited state are degenerated. However, this assumption was found 
to work well in a more precise calculation with the spin of t and the spin-dependent 
a-t interaction taken into account explicitly. Namely, the S'-factors of the transition 
to the ground and excited states were found to be, respectively, 53% and 50% of 
the ^-factor by the present spin-neglected calculation (the solid curve in Fig. 3) in 
the Gamow peak region, and therefore their sum (103%) deviates from the spin- 
neglected case by only 3%, although the above respective percentages deviate from 
the numbers 67% and 33% (proportional to the spin weights) in the case where the 
3/2" ground and l/2~ excited states are assumed to have the same degenerated 
binding energy and the same radial wave function. 

The magnitude of S{E) for ^Li production in CBBN is approximately 30 times 
larger than the S*- factor of the El radiative capture SBBN reaction (2-6) in the 
Gamow peak region in Fig. 3, Sj{Eo) ~ 0.1 keV b.^^) This CBBN/SBBN enhance- 
ment factor of ~ 30 is very much smaller than that of ~ 10^ for ®Li production. 
This is simply because the El transition in SBBN is allowed in ^Li production but is 
heavily hindered in ^Li production. However, as was pointed out in §2.1.2., Cyburt 
et al.^^ predicted a very large CBBN/SBBN enhancement factor (~ 10^) assum- 
ing the scaling relation.^) This large overcstimation is because the scaling relation, 
which compares the wavelengths of the real and virtual photons, is not suitable for 
this type of transfer reaction caused by the strong nuclear interaction. On the other 
hand, Kusakabc et al}"^^ neglected the CBBN processes for the creation of ^Li and 
'^Bc in their BBN network calculation by taking the same consideration as above for 
the angular-momentum rearrangement (^i = Li = 0) — > (I2 = L2 = 1). Although 
consideration of the largest contribution from {li = 0, Li = 1) ^ (^2 = 1,-^2 = 0) is 
missing in Rcf. 12), we support their assumption as a reasonable one. However, our 
reaction rates given below will be employed in the BBN network calculation when 
the abundances of ^Li and ^Be are precisely discussed. 

The energy dependence of S{E) in Fig. 4 may be approximated by expression 
(2-27) with S{0) = 2.6 keV b and a = 0.02 b, and therefore the reaction rate is 
written, for Tg < 0.2, as 



NA{av) = lAx 10'^ Tg 3exp(-6.08r9 + l.STg^ + 0.55Tg) cm^s'^mor^ 

(2-38) 

The calculated S'-factor for the ^Be production reaction (2-5), is shown in Fig. 4. 

Almost the same discussion as that for ^Li production can be made. The 5-factor 
S{E) in Fig. 4 may be approximated with S{0) = 13.7 keV b and a = 0.01 b, then 
the reaction rate is written, for Tg ^ 0.2, as 



2 12 

Na{(tv) =9.4x 10^T9~3exp(-9.66r~3)(l + 0.20r93 +0.05r9) cm^s^^mor^ 

(2-39) 
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Fig. 4. The 5- factor of the CBBN reaction (2-5) obtained by the three-body calculation (solid line). 
The S-, p- and d-wave contributions are shown individually. The dotted curve (in arbitrary units) 
illustrates the Gamow peak for Tg = 0.1 (fcT = 8.6 keV) with Eo = 60 keV. 



§3. X -catalyzed radiative capture reactions 

3.1. Necessity of three-body calculation 

Since the A = 7 nuclei dominantly produced are ''Be (eventually decaying to 
^Li by electron capture) in the BBN network calculation with the CMB-based rjB = 
6.0 X 10~^°, any reaction*) that destroys ^Be might be effective in reducing the 
overproduction of ^Li. For this purpose, Bird et al.^^ considered a resonant radiative 
capture reaction (cf. Fig. 5), 

CBeX-) {^BX-)^^f ^ (8bX-) + 7(~ 0.7 MeV) , (3-1) 

where the intermediate state, denoted as {^BX~)2p^\ is a Feshbach resonance gen- 
erated by the coupling of the atomic 2p excited state, say {^BX'~)2p, with the 
('''BeX~)-|-p continuum. Here, note that after rapid /3-decay, the bound state (^BX~) 
transforms to (^Be(2+, 3MeV)X~), which immediately decays**'' to the scattering 
channel a-\-a + X' + 1.5 MeV. 

Assuming a simple Gaussian charge distribution of the ^B nucleus. Bird et al. 
calculated the energy of {^BX~)2p, E2p, with respect to the ^B+X" threshold. They 
obtained E2p = —1.026 MeV and estimated the resonance energy, £^res) with respect 
to the ('^BeX") +p threshold as E^^s = E2p + Es-b — Ef^r^ex-) = 0.167 MeV without 

The reaction CBeX-)+p -^^B + X' is energetically impossible due to the negative Q- value 
(-1.323 MeV); note that the binding energy of *B is only 0.138 MeV. 

**•* The authors of Rofs. 6) and 14) erroneously consider that (*BX~) transforms to (*Bc(gs)X~) 
and that the latter state could potentially lead to a new primordial source of '""Be via (**BeX^) + n 
-> ^Be + X". Note that ^B shows no /3-decay to ^Be(gs). 
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Fig. 5. Schematic illustration of the X -catalyzed resonant and nonresonant radiative capture 
processes (3-1) and (3-2). 



the scattering calculation. Further assuming that the resonance width i^es is much 
larger than the radiative width for decay to {^BX~)is, they derived the reaction 
rate of (3-1) and claimed that the resonant reaction (3-1) is effective in reducing the 
^Li-^Be abundance in their BBN network calculation^-* and in restricting the lifetime 
and primordial abundance of the X~ particle. 

As discussed below, the rate of the resonant reaction depends strongly on Ej-cs, 
but the model in Ref. 6) is too simple to treat the resonance state. Therefore, in 
this section, we perform a precise ^Be + p + X~ three-body scattering calculation of 
the resonant reaction (3-1) as well as the nonresonant radiative capture reaction 



(^BeX-)+p^ (SBX-)+7(> 0.6MeV), 



(3-2) 



although the reaction rate of (3-2) is much smaller than that of (3-1). 

Since the rate of the resonant reaction is proportional to exp(—E,.(^/kT), a small 
change in E^es can generate a large change in the rate; for example, at Tg = 0.3 {kT = 
25 kcV), increase (decrease) of 50 keV in Ej-^s changes the reaction rate by a factor of 
e^(= 7.4). The 50 keV change is only 5% of the binding energy {—E2p ~ 1.0 MeV) 
of the {^BX~)2p state, and therefore it is essential to calculate Ej-^s with the detailed 
structure of the ^B nucleus taken into account. 

Since the nucleus ^B has been extensively studied, mainly from the viewpoint 
of the physics of unstable nuclei and that of the solar neutrino problem, a large 
amount of information about ®B has been accumulated.^^)' The ground state 
of ^B, being very weakly bound (^sb = —0.138 MeV) with respect to the ^Be + p 
threshold, is known to have a p- wave proton halo (a long-range tail) with the '''Be-l-p 
structure. Therefore, as will be shown below, the charge distribution of ^B has 
a long-range quadrupole component, which generates the anisotropic part of the 
Coulomb potential between ^Be and X~. The expectation values of this part for 
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the {^BX~)2p wave function amount to some —50 keV (J = 0), +10 keV (J = 1) 
and +30 keV (J = 2), where J denotes the total angular momentum, to which the 
p-wave ^Be — p relative motion in and the p-wsve relative motion between and 
are coupled (see (3-3) below). 

Another important factor that we have to consider concerning ®B is the spins 
of the ^Be (3/2~) core and the valence proton (1/2). The two spins couple with the 
p-wave angular momentum between ^Be and p to give the total angular momentum 
of 2"*" in the ground state. The other possible spin-coupling states are not bound; 
the excitation energies of the 1^" and 3"*" resonance states are 0.77 and 2.32 MeV, 
respectively. Therefore, we cannot neglect the fact that the ground state of ^B has 
the specific spin of 2+. 

However, the estimation of E2p in Rcf. 6) was made assuming an isotropic Gaus- 
sian form of the charge distribution of *^B with neither the ^Be+p structure nor the 
spins of the particles taken into account. In this section, we investigate reactions 
(3-1) and (3-2), explicitly adopting the '''Be + p + X~ three-body degree of freedom 
in which the very diffuse, anisotropic charge distribution of ^B (= ^Be + p) is au- 
tomatically taken into account. The use of this three-body model makes it possible 
to calculate not only the precise resonance energy £^res but also the cross section 
of the first transition process in (3-1), which could not be treated by the model of 
structureless ^B in Ref. 6). 

3.2. Resonance energy obtained using approximate models o/^B 

Before performing the three-body calculation, we employ four types of approx- 
imate models of ®B, Models i) to iv) below, and discuss how the energy of the 
{^BX~)2p state, E2p, depends on the assumed structure of ^B, although we do not 
couple the state to the ^Bc + p scattering state. This may be an instructive guide to 
the sophisticated three-body scattering calculation in §§3.3 and 3.4. 

Model i) Bird et alS'^ assumed a Gaussian charge distribution of ^B with an 
rms charge radius of 2.64 fm. This gives i?2p = —1.026 MeV {mx ^ oo is taken 
throughout this subsection) , which corresponds to a resonance energy of = 167 
keV above the (^BeX~)+p threshold at —1.330 MeV with respect to the ^Be+X^+p 
three-body breakup threshold (see Fig. 5). however, this assumption of the Gaussian 
charge density is not appropriate for this special nucleus ^B. 

Model ii) The ®B nucleus is known to have a very loosely bound ^Be+p structure 
with a ^>-wave proton halo around the ^Be core.^^^"^^) Neglecting the spins of ^Be 
(3/2~) and the valence proton (l/2~), we firstly calculate the p-state ^Be-p wave 
function using the most well known ''Be-p potential^^) (parameters are given in §3.3). 
By assuming the charge density of ''Be (proton) to be the Gaussian with the observed 
rms charge radius of 2.52 fm^^^ (0.8750 fm^^^), we calculate the charge density of 
the ground state of ^B. Owing to the p-state wave function, the density has both an 
isotropic monopole part and a deformed quadrupole part. The former is illustrated 
in Fig. 6 together with the contributions from ^Be and p. The density significantly 
deviates from the Gaussian shape in the tail region owing to the proton halo. If 
we take the ^B-X~ Coulomb potential due to the monopole charge density alone. 
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we have £;(™°°°)(2p) = -1.009 MeV, which corresponds to E^} = 184 keV. This 
increase of E^es by 17 keV, compared with 167 keV in Model i), reduces the reaction 
rate of (3-1) by a factor of 2 at Tg = 0.3. 




r[fm] 

Fig. 6. Monopole part of the charge density, p{r), of (dotted line) caleulated using the '^Be + p 
model. The solid (dashed) line shows the contribution from the proton (^Be). p(r) is normalized 
as / p{r)dr = 5, r being the distance measured from the cm of *B. 



Model iii) Using the above Model ii), wc further consider the contribution from 
the quadrupole part of the charge density of ^B. Let 4)\m{f2) and iP\m{R2) denote 
the wave functions of the p-wave ®B ground state and the atomic 2p relative motion 
between and X~ , which have already been obtained in Model ii). Then, 

•^^S = ['^1(^2) ® MR2)]jM {J = 0, 1,2) (3-3) 

is the total wave function of the (^BX^)2p state*^ with angular momentum J and 
z-component M. The quadrupole part of the Coulomb potential between *B and 
X~ , which is sensitive to the angle between r2 and -R2, should contribute to the 
total energy dependently on the angular momentum coupling to J. The expectation 
value of the Hamiltonian with respect to ip'j^ , say -Bj"^ (2p) , is written as 

Ef^ {2p) = {2p) + AEf'^'^ {2p), (3-4) 

where /\E^^^^'^^ (2p) is the contribution from the quadrupole part of the charge den- 
sity of ^B mentioned above and (2p) was given in Model ii). Wc obtain 
^^Cquad)^2p) = -54keV (J = 0), +27keV (J = 1) and -5.4 keV (J = 2), and there- 
fore, Ef'\2p) = -1.063keV(J = 0), -1.036keV(J = 1) and -1.014keV(J = 2), 

*^ In this notation of (*BX~)2p in Model iii), and similarly in Model iv), the p-wave angular 
momentum of *B and the total angular momentum J are not explicitly written for simplicity. 
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which correspond to E^r^es^ = 130 keV (J = 0), 211 keV (J = 1) and 179 keV (J = 2). 
In the case of J = 0, the resonance energy, Ey-^s = 130 MeV, results in a reaction 
rate that is 8 times larger than that of Model ii). However, we have not yet included 
the effect of the spins of ^Be (3/2-), p (1/2) and (2+). 

Model iv) The spin-parity of the ground state of ^B is 2^. In previous studies on 
^B taking the ''Be + p model, the structure of low-lying states of ^B is considered as 
follows: The valence proton is located in the orbit around the ''Be (3/2^) core, 
and the two 3/2" spins are coupled to the total angular momentum / = 2+ in the 
ground state and / = !"'" and 3^ in the excited states with the observed = 0.77 
and 2.32 MeV, respectively. Since the spin-dependent interaction is strong and the 
Coulomb interaction between ^B and X~ does not change the spin structure, we 
have to seriously consider the ground-state spin 2+ and its coupling scheme in the 
calculation of the resonance energy. Let 4>{p.^^,^^s/2}iMi{''^2) denote the wave function 
of the ground state (/ = 2) of ^B explained above. The total wave function of the 
{^BX~)2p state is written as 



</'{p3/2,3/2}/(r2) ® MR2)\ , (3-5) 



where V'i(-'^2) was given in Model ii). The expectation value of the Hamiltonian 
with respect to ^j]}, say E^J^j\2p), is described by the form 

2 

4})(2p) = E(-°-)(2p) + J2 aj,iAAE^r^\2p), (3-6) 

where E^'^°'^°\2p) and AE^^^'^\2p) are given in (3-4), and ajjA is written as 

1 ^ 
aj,iA = ^{2J + 1){2A + 1) ^ [W{IIAJ- Sl)f . (3-7) 

The second term of Eq. (3-6) vanishes owing to the summation over S for the ground- 
state spin / = 2 but not for I ^ 2. We then have Efj^{2p) = E^'^'"^°\2p) = -1.009 

MeV for / = 2+, and therefore E^Z' = 184 MeV (J = 1, 2) which is, by chance, the 
same as in Model ii). 

3.3. Three-body calculation of resonance energy and width 

In the four types of approximate models for {^BX~)2p in the previous subsection, 
^B is assumed to be the same as that in the free space. However, since ^B is a very 
loosely bound state of ^Be and p, it is possible that the spatial structure of ^B 
changes in the presence of X^, namely, in the presence of the Coulomb potential 
between ^Be and X~ and that between p and X~ (the spin structure is not changed 
by these Coulomb potentials). Therefore, using the ^Be + p + X~ three-body model 
with the Hamiltonian (2-9), we calculate the resonance state as a Feshbach resonance 
embedded in the (^BeX~) +p continuum and precisely determine the energy and 
width of the resonance. 
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3.3.1. Nuclear and Coulomb potentials 

In the ^Be + p model, we treat ^Be as an inert core with spin 3/2~. We assume 
Gaussian charge distributions of ^Be and the proton as 4e(7r6|g)~^/^e~(^/''Be) and 
e(7^6p)~^/^e~('■/^^')^ respectively, and take 6Be = 2.06 fm and bp = 0.714 fm to 
reproduce the observed rms charge radii, 2.52 fm^^^ for ^Be and 0.8750 fm^^^ for the 
proton. The Coulomb potential between ^Be and X~ is then given by 

V,^,ir) = -Ae''EMM, (3.8) 



and that between X and p is written as 

Vp-x{r) = -e' (3.9) 



2 erf (r/6p J 



The energy of {^BeX-)is is = -1-386 (-1.324) MeV and the rms radius is 3.60 
(3.49) fm for mx = 100 GcV (mx 00). 

The potential V7-Bc-p{f) is a sum of the nuclear potential, VT^^pir), and the 
Coulomb potential, V'f^^^{r). The latter is given by 



erf(r/^6L + 6|) 



^r^e-pW = ^ -. (3-10) 

For the nuclear potential ^Be-p(^) between ^Be and p, we follow the work of Ref. 38), 
which is a standard study on the ^Be + p ^ ^B + 7 reaction. The nuclear potential 
is parameterized as a Woods-Saxon potential plus a spin-orbit potential with an 
adjustable depth Vb(Zj, /): 



Fso{l • s) — — 
r dr 



l + exp[(r-i?o)/a] 



Here, we take a = 0.52 fm, Rq = 2.391 fm and Fso = 0.351 fm.^^) The orbital 
and total angular momenta of the proton are denoted by (Ij). The 1 = 2+ ground 
state of ^B is described in terms of a pure (Ij) = P3/2 orbit of the proton coupled 
with the 3/2" spin of the '''Be core. The well depth Vq {I = l,j = 3/2,1 = 2) is 
adjusted to reproduce the binding energy of 0.1375 McV of the ground state and 
is given by —44.147 MeV. The s-wave potential between ''Be and p is assumed to 
have Vq = —30.0 MeV, Rq = 2.39 fm and a = 0.65 fm and additionally have a Pauli 
repulsive potential with Vq = 800 MeV, Ro = 1.5 fm and a = 0.2 fm. This s-wave 
potential gives no bound state. 

3.3.2. Resonance wave function 

Let ^3 (^Bc) and ^1 (p) denote the spin functions of the ''Be ground state and 

the proton, respectively. The Is ground state of (''BeX") is described by ^qq (^"1). 
We consider the s-state relative wave function between the incident proton and the 
target (^BeX~), say Xto (-Ri)(= Xo\Ri)^oo{Ri)) , since that of the p-state has a 



Big-Bang Nucleosynthesis Reactions 



21 



different parity and that of the d-state has neghgible contribution to the resonant 
radiative capture process. There are neither inelastically excited channels in c = 1 
nor transfer channels in c = 2 and 3 for the energies (< 200 keV) concerned here. 
The total wave function used to describe the resonant state is then written as 



<^^jM = tfJ(n)xS(i2i) ^3rBe)®ei(p) „ +^yM"'^- (3-12) 



^(closed) 



JM 



The second term of (3-12), ^^^^^'^\ represents the internal amplitude of the 
resonance wave function, whose approximate expression is ^j^j in (3-5). However, 

we describe ^j^^^'"^'' taking the three-body degrees of freedom in the same way as in 
§2.4. Since we find that, in (3-12), the contributions from the amplitudes with c = 1 
and 3 are negligible in the resonance energy region, we here express ^j^^^^^^ as 



^(closed) ^^^^^^(2)^^^ (3-13) 



.(closed) _ ^ J, - ^(2) 
i/=l 

using the three-body basis functions only in c = 2 together with the spin functions, 

^jL,u =T.^%n.,N, [[ [^n,l,{r2) ® ® ® (^^2)] , (3-14) 

where I2 = L2 = 1,1 = 2 and J = 1,2 are sufficient to describe the resonant 
reaction. For the Gaussian ranges in (2-34) and (2-35), we take nmax = -^max = 15 
(^'max = 225) and {ri, r„,,^^^^, i^Ar^.^,} = {0.4, 15.0, 0.6, 20 fm}, which is sufficiently 
precise for the present purpose. 

By diagonalizing the Hamiltonian, we obtain eigenstates {^jMy, ly = I — t'max}, 
among which the energy of the lowest-lying state {^jm,u=i) measured from the 
'^Be+p threshold is 198, 186, 177 and 174 kcV for mx = 50, 100 and 500 GeV and 
mx — 00, respectively (J = 1,2). The mx dependence of the energies comes from 
the fact that the kinetic energy and the "^Be + p threshold energy depend on mx- 

The total wave function (3-12) is solved under the scattering boundary condition 

lim RiX^o\Ri) = ui-\h,Ri) - S(^^ui+\h, Ri) (3-15) 

itl— >0O 

on the basis of the same prescription as that in §2.2. The resonance state should 
appear around the energy of the pseudostate ^jm, v=i mentioned above. The calcu- 
lated partial- wave elastic scattering cross section is illustrated in Fig. 7 for J = 1~ 
with mx = 50, 100 and 500 GeV and nix oc; similar behaviour is obtained for the 
resonance with J = 2~ . The energy and the proton width of the resonance 
arc summarized in Table I. 

3.4. Result for the resonant radiative capture 

The final state, {^BX^), of the radiative reactions (3-1) and (3-2) is obtained as 
the ground state of the '^Be+p + X~ system. The dominant component is obviously 
the product of the 2+ ground-state wave function of and the Is wave function of 
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Fig. 7. Partial- wave (J = 1 , s-wave) cross section of the elastic CBeX)is +P scattering near 
CBeX-y^"; for the cases of mx = 50 GeV (dash-dotted hne), 100 GeV (solid line), 500 GeV 
(dashed line) and mx — > oo (short dashed line). Similar behaviour is seen for J = 2~ (s-wave). 



(^BX^). Wc describe it more precisely using the same three-body basis functions in 
(3-14) with I2 = 1, L2 = 0, 1 = 2 and J = 2 with the same ranges of the Gaussian 
basis. Other configurations are not necessary within the accuracy required for the 
present purpose. By diagonahzing the three-body Hamiltonian, we obtained the 
ground-state wave function, ^j^2+M energy Egs, which is hsted in the seond 

to last column of Table I. 

The width of the electric dipole (El) transition from the resonance state ^j^^ 
with J~ to the 2+ ground state is given by 

^' = ^2jTT ^ K«'I<^I•^^M^)I^ (^ = 1'2) (3-16) 

MM' n 

where the El operator Qi^^^ is defined by 

QSf)= J^g,i?pyi^(^f). (3-17) 

i=l 

Here qi and Rf are the i-th particle's charge and its position vector relative to the cm 
of the total system, respectively, and A;^ is the photon wave number. The calculated 
values of 7"!^ ^^'c listed in Table I; they are, as expected, close to the width (10.4 eV) 
given by the simple atomic El transition {^BX~)2p {^BX~)is- 

Finally, the reaction rate of the resonant radiative capture process is given using 
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Table I. The '^Be + X~ + p three-body calculation of the energy (Eres), the proton decay width 
(Fp) and the radiative decay width (F-y) of the resonance states with J = 1~ and J = 2~ as well 
as the energy (-Egg) of the three-body ground state with J = 2^. Eres and -Egs are measured 
from the C'BeX~)is + p threshold, whose energy (Sth) is given in the last column with respect 
to the three-body breakup threshold. All results are calculated for mx = 50, 100 and 500 GeV 
and mx — > oo. 



mx 
[GeV] 




J =1- 


(res) 




J = 2- 


(res) 


J = 2+ (gs) 

Egs 
[keV] 


threshold 

Eti, 
[keV] 


E'' 

-'-'res 

[keV] 


[keV] 


[eV] 


E-' 

-'-'res 

[keV] 


[keV] 


[eV] 


50 


196.4 


0.90 


9.1 


196.9 


0.55 


9.1 


-624.2 


(-1252.0) 


100 


185.0 


0.82 


9.6 


185.5 


0.50 


9.6 


-635.5 


(-1286.1) 


500 


175.8 


0.74 


9.9 


176.3 


0.44 


9.9 


-643.6 


(-1316.4) 


oo 


173.0 


0.71 


10.1 


173.6 


0.43 


10.1 


-645.9 


(-1324.0) 



the formula for the reaction (Eq. (4-194) in Ref. 26)) as 

^-^^^) = ^-^ (m^J g (2/, + i)(2i. + i)r/+r,^ -H"^J' 

(3-18) 

where Ii and I2 arc the spins of ''Be (3/2) andp (1/2), respectively. Here, we consider 
that the total width of the resonance, i^es) is given by Fp + since there is no 
other decaying channel. We replace the resonance energies E^^{J = 1,2) by their 
average, since they are almost the same. We then obtain (in units of cm^ s~^ moP^) 



NA{av) 


= 1.37 X 10^ Tg' 


^ exp (- 


-2.28/T9), 


{mx = 


50 GeV) 


(3-19) 


Na{(7v) 


= 1.44 X 10^ Tg" 


.3 

^ exp (- 


-2.i5/r9), 


{mx = 


100 GeV) 


(3-20) 


Na{(7v) 


= 1.48 X 10^ Tg" 


3 

^ exp (- 


-2.04/r9), 


{mx = 


500 GeV) 


(3-21) 


Na{c7v) 


= 1.51 X 10^ Tg" 


3 

^ exp (- 


-2.oi/r9). 


{mx 


00) 


(3-22) 



The mx dependence of the rates is not negligible since, at Tg = 0.3, their ratio 
is 0.37 : 0.60 : 0.88 : 1.0 for mx = 50, 100 and 500 GeV and mx ^ 00. The rate 
given by Bird et al. in the second part of Eq. (3-22) in Ref. 6) for mx 00 is by 
chance close to the above rate (3-22); at Tg = 0.3, their rate is 1.2 times that of ours, 
although the models are markedly different from each other. 

3.5. Result for the nonresonant radiative capture 

In this subsection, we investigate the nonresonant (direct) radiative capture 
process (3-2) using the '^Be+p + X~ three-body model. Here, we ignore the intrinsic 
spins of ^Be and p because, as will be shown below, the calculated reaction rate 
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without the spins is 4 orders of magnitude smaher than that of the resonant reaction 
(3-1). The further inclusion of the spins will not change the result meaningfully. 

We consider the nonresonant El radiative capture from the J = incoming 
state of the (^BeX^) + p channel to the J = 1 ground state in which the ^Bc and p 
are dominantly in p-wme relative motion (other partial- wave states are negligible), 
and we denote the wave functions respectively as l^oo(-^^) and ^^^^ with the proper 
normalization. The wave functions can be obtained using the same methods as those 
in the previous subsections, but the second term of (3-12) can be neglected for the 
present nonresonant scattering wave. 

The cross section of the El capture is given by 



.(El) 
cap 



(E) 



167r fe3 



(3-23) 



M 



where is the wave number of the emitted photon and vi is the velocity of the 
relative motion between (^BeX~) and p. 

The calculated S{E) of the CBBN reaction (3-2) is illustrated in Fig. 8 together 
with the observed S-factor^^^ of the SBBN partner reaction 



'^Be+p^ ^B + 7(> O.UMeV). 



(3-24) 
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Fig. 8. The astrophysical S-factor (solid line) of the nonresonant radiative capture CBBN reaction 
(3-2) obtained by the three-body calculation. The observed S-factor of the SBBN partner 
reaction (3-24) is shown by open circles.^*' The dotted curve illustrates the Gamow peak (in 
arbitrary units) for Tg = 0.3 (kT = 35 keV) with the maximum at Eo = 114 keV. 



The enhancement ratio CBBN/SBBN of the S-factor is only 1 - 2 in the Gamow 
peak region for the nonresonant radiative capture processes (3-2) and (3-24). How- 
ever, in the work of Bird et al.f^ this ratio is estimated to be as large as ^ 700. 
The origin of this enormous overestimation is the crude assumption employed in 
their model. Namely, they assumed that the magnitude of the three-body El matrix 
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element in Eq. (3-23) is the same as that of the two-body El matrix element in 
SBBN. This is an incorrect assumption because, in the SBBN matrix element, the 
dominant contribution comes from the asymptotic tail region of the distance between 
^Be and the proton (along in the loosely bound ground state of ^B, but, in the 
CBBN matrix element, the contribution is heavily suppressed by the presence of X~ 
as follows: (i) in the (^BeX~)is state of the initial channel, the distance between 
'''Bc and along ri is short (the rms radius is 3.5 fm), (ii) in the final state 
!P'oo(£'), the distance between ^B and X^ along R2 is also short (the rms radius is ~ 
3 fm), and thus (iii) these strongly confine the possible ^Be-p distance along r2 that 
is effective in the three-body El matrix element. Therefore, the contribution from 
the asymptotic region along r2 is heavily suppressed. The CBBN/SBBN ratio of 
the S'-factor becomes larger in the scaling model of Rcf. 6) as the binding energy of 
(f'BX^) increases, but, at the same time, the wave function of C^BX^) shrinks and 
therefore the three-body El matrix element becomes smaller. This consideration 
of this mechanism was missing in the simple model of Ref. 6) for the nonresonant 
radiative capture reaction. 

The fact that the CBBN and SBBN S'-factors have similar magnitudes, within 

factor nearly 2, in the Gamow peak region in Fig. 8 suggests an approximate model 

/pi \ 

for the cross section Ucap {E) of the CBBN reaction. The cross section may roughly 
be derived from (2-21) by employing the observed S-factor of the SBBN-reaction 
and the Coulomb barrier penetration factor of the CBBN channel. Here, the SBBN 
partner reaction is defined as the reaction that is given by simply removing X" from 
a CBBN reaction*^ This approximate model will be further examined in the next 
section for complicated CBBN three-body breakup reactions. 

In Fig. 8, since the energy dependence of S{E) may be approximated by (2-27) 
with S{0) = 20 X 10~^ keV b and a = 0.15 x 10~^ b, the reaction rate is written as 

NA{av) =2.3x 10^T9"^exp(-8.83r9"^)(l-M.9r/ -F0.54r9) cm^s-^mol-\ 

(3-25) 

for Tg < 0.5. Since this rate is 4 orders of magnitude smaller than that of the 
resonant reaction (3-1) at Tg = 0.3 — 0.5, the nonresonant radiative capture reaction 
(3-2) will play a very minor role in the BBN network calculation. 

3.6. Comment on resonant recombination between ^Be and X~ 

The radiative capture processes (3-1) and (3-2) are preceded by the recombina- 
tion of ^Be and X~ to form (^BeX~). Bird et al.^^ calculated the reaction rates 
of the resonant and nonresonant processes of the recombination. We here make an 
important comment on a part of their calculation. In addition to the normal re- 
combination processes (2-7) and (2-8) in Ref. 6), they specially considered another 
resonant recombination process, (2-12) in Ref. 6), 

^Bes- +X-^ CBe,-X-)2s ^ C^e,-X-)2p + 7 ^ (^Be3-X-)i, + 37, (3-26) 
222 2 

*^ For instance, the reaction a + d — > ®Li + 7 is not the SBBN partner of {aX~) + d ^ ®Li + X~ 
but that of (aX") + d ^ {^UX') + 7. There is no SBBNpartner of (aX") ^Li -|- X" . 
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where Be 3- and 'Bei- are the ground state and first excited state at = 0.429 

2 2 

MeV, respectively. Their calculated energy of the intermediate state (^Bei-X^)2s 

2 

with nix ^ oo is 49 keV below the ^Bca- + threshold, which means that the 

2 

intermediate state cannot be a resonance in (3-26). Bird et al., however, artificially 
pushed the state upward to a resonant energy of +10 keV above the threshold and 
derived a large recombination rate, Eq. (2-13) of Ref. 6). Their reason for doing 
this was that they expected 'a sizable nuclear uncertainty' such as a correction for 
the finite mass of , a larger charge radius of the excited state and a correction 
for nuclear polarizability. The contribution of the 'resonant' (^Bei-X~)2s state is 

2 

then capable of enhancing the recombination rate by a factor of a few. Bird et al. 
reported two types of results, with and without the (^Bei-X~)2s contribution, when 

2 

they discussed the abundance of ^Li-^Be and the lifetime and abundance of X~ 
particles. 

To completely remove the above 'nuclear uncertainty' from the model, we per- 
form an a + ^He + X" three-body calculation of the energy of (^Bei-X~)2s. The 

2 _ 

three-body Hamiltonian is the same as that used in §2, where we investigated the re- 
action (aX^)+'^He ^Be + X~, but here we take into account the spins of ^He and 

^Be using a spin dependent ct-^He potential to reproduce the energies of ''Be 3- and 

2 

'''Be 1 - . We find that the rms charge radius differs by only 0.05 fm between '''Be 3 - and 

2 2 

^Bei-. According to the three-body calculation the energy of (^Bei-X^)2s is —41 

keV^for mx ^ 00, -20 keV for mx = 100 GeV and -2 keV for mx'= 50 GeV with 

respect to the ^Bea- +X~ threshold. Therefore, we conclude that the (^Bei-X~)2s 

2 2 
state never becomes a resonance in the process of (3-26). In Ref. 6), two types of 

calculations of the element abundance with and without the (^Bei-X~)2s resonance 

2 

are reported, but the case with the resonance is not acceptable. 

The application of the a + ^He + X^ three-body calculation to the full recombi- 
nation processes of ^Be and X" will be one of our future subjects of study since the 

transition between the ^Bca- and ^Bei- states due to X~ is an important factor in 

2 2 

the recombination and can be unambiguously treated using the three-body model. 
§4. X~-catalyzed three-body breakup reactions 

The most effective reactions for destroying ^Li and ^Li in SBBN are 

^Li + p^ a + =^He + 4.02MeV, (4-1) 
^Li + p^ a + a + 17.35 MeV, (4-2) 

which have large S'-factors (5(0) ~ 3 MeV b and ~ 0.06 MeV b, respectively^^)). 
Therefore, the corresponding CBBN three-body breakup reactions, 

{^UX-)+p^ a + ^He + X- + 3.22MeV, (4-3) 
CuX-)+p^ a + a + X- + 16.47MeV , (4-4) 
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should be taken into account in BBN calculations* ^ . However, the explicit calculation 
of these CBBN processes is tedious and difficult because (i) the exit channel is 
that of a three-body breakup, (ii) at least a four-body model, in which ^Li ('^Li) is 
composed of d + a {t + a), is needed, (iii) d{t) in ^Li (^Li) must be picked up by the 
incoming proton to form ^He (a) and (iv) it is tedious to reasonably determine all the 
three spin-dependent nuclear interactions appearing in the four-body model. Instead, 
a naive approximation often taken in BBN network calculations, for instance, in 
Rcfs. 4) and 12), is that the cross section of the CBBN reaction (4-3), and similarly 
that of (4-4), may be given by a product of the observed 5- factor of the SBBN 
partner reaction (4-1) and the CBBN Coulomb barrier penetration factor, in which 
the charge of the target is reduced by one unit owing to the presence of X~. 

In this section, we propose a more sophisticated and phenomenologically reason- 
able three-body model to implement the information of the SBBN cross section into 
the CBBN calculation instead of carrying out an explicit calculation of (4-3) and 
(4-4). Firstly, we calculate the SBBN reactions (4-1) and (4-2). We do not explicitly 
treat the channel coupling between the entrance and exit channels. Instead, we em- 
ploy only the entrance channel, say {AX~) + a, and introduce a complex potential 
VA-a{f) between the particles A and a (here, A =^'^Li and a = p): 

VA-air) = vlf{r) + iV;i-^\r) (4-5) 

as seen in nuclear optical model potentials. In this case, the absorption cross section 
in the elastic A+a scattering is equivalent to the reaction cross section because there 
are no open channel other than the entrance and exit channels of (4-1) and (4-2). We 
determine the potential VA-a(r) so as to reproduce the observed SBBN cross section 
(5-factor). 

Secondly, we incorporate this potential VA-air) into the three-body Hamiltonian 

(2-9) of the A + a + X~ system and solve the elastic scattering between {AX~) 
and a, namely the elastic (®''^LiX~) +p scattering. We consider that the absorption 
cross section obtained in this scattering calculation provides the cross section of the 
CBBN reaction**). 

The wave function is written similarly to (2-12) but without the exit channel as 

^JM = 0S (r-i) x%{Ri) + , (4-6) 

where ^^oo^(''i) represents the Is wave function of (®'^LiX~) and Xjm(-Ri) for the 
(^'^LiX^) 4- p scattering wave. The scattering boundary condition imposed on 

X?i(iii)(= X^PiRi)YjMiRi)) is given by 

hm Rix[i\Ri) = u[j-\ki,Ri) - S(^^u[j^\h,Ri). (4-7) 

Ki— »oo 



Since the kinetic energy of the exit channel is large, the formation of the bound states '^Be, 
®Be, {aX~) and (^HeX~) in (4-3) and (4-4) is not important and is not explicitly considered there. 

**•* More precisely, this absorption cross section, for instance, in the case of {^LiX~)+p scattering, 
includes transitions to the {aX~) + ^Hc, (''HeX") + a and '^Be + X~ channels, which are possible 
in the presence of the X^ particle, but are of minor importance owing to the much smaller phase 
space than that in the three-body breakup channel. 
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Similarly to in the previous sections, the second term of (4-6), •f'j^**^'^'* , represents 
all the asymptotically vanishing three-body amplitudes that are not included in 
the first scattering term, and is expanded in terms of the eigenfunctions of the 
Hamiltonian (without the imaginary part of (4-5)) as 

^!jT''^ = Ebj-^^jM,u. (4-8) 
i/=i 

Using S(^-^ in (4-7), we derive the reaction cross section as 

oo 
'^1 J=0 

This reaction (absorption) cross section can be expressed alternatively as 

^Treac = ^^{"^JM I^i^f''^(r2) Wjm)- (4-10) 

These two types of o"rcac utilize information from different parts of the three-body 
wave function, namely, the information from the asymptotic part along Ri in the 
former expression and that from the internal part along r2 in the latter. Therefore, 
it is difficult to demonstrate the agreement between the two types of (Tj-eac- 

We 

obtained a precise agreement between their numbers to four significant figures, which 
demonstrates the high accuracy of our three-body calculation. 

4.1. Nuclear complex potential and Coulomb potential 

We construct the complex potential*^ (4-5) between ^'^Li and p so as to reproduce 
the observed low-energy 5-factors of (4-1) and (4-2) (the adopted values in Ref. 18)). 
The shape of the potential is assumed as 

y(reai)^^) = Vo [1 + exp{(r - Ro)/ao}]-^ + VLip(r), (4-11) 
y(imag)(^) = Wo [1 + exp{(r - Rj)/ai}r^ . (4-12) 

The Coulomb potential V^^_p{r) is derived by assuming a Gaussian charge distribu- 
tion for ^Li (^Li) that reproduces the observed rms charge radius 2.54 fm (2.43 fm).^^) 
The energies of {^LiX~) and (^LiX~) are -0.773 and -0.854 MeV, respectively, and 
the rms radii are 4.61 and 4.16 fm, respectively, when mx = 100 GeV. 

i) ^Li-p potential 

The potential parameters are chosen as follows: Vq = —59.5 MeV, Rq = 2.5 fm 
and a = 0.6 fm for the odd state and Vq = 400 MeV, Rq = 1.5 fm and a = 0.2 fm 
for the even state as well as Wq = —5.0 MeV, Rj = 2.5 fm and aj = 0.6 fm for 
both states. The repulsive potential for the even state is introduced as a substitute 
for the Pauli principle to exclude the s-state (the contribution from the d-state is 
negligible) . 



*^ We follow the same prescription as that employed in Ref. 40) to determine the nuclear fusion 
rate in a muonic molecule {dtfj,) by using a complex potential between d and t. 
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ii) ^Li-p potential 

The potential parameters are chosen as follows: Vq = —77.2 MeV, Rq = 2.5 fm 
and a = 0.6 fm for the odd state, Vq = 400 MeV, Rq = 1.5 fm and a = 0.2 fm for 
the even state, and Wq = —1.1 MeV, Ri = 2.5 fm and a/ = 0.6 fm for the odd state 
and Wo = for the even state. The reason for setting Wo = for the even state is 
that, in reaction (4-4), the odd angular momentum between the two a is prohibited, 
and therefore the even state of the ^Li-p relative motion is forbidden ('''Li has odd 
parity). This explains why the observed S'-factor of (4-4), <S'(0) ^ 0.06 MeV b, is 
much smaller than that of (4-3), S{0) ^ 3 MeV b. 

4.2. Result for three-body breakup reactions 

The calculated S'-factors of the CBBN reactions (4-1) and (4-2) are shown in 
Figs. 9 and 10, respectively, together with those of the SBBN partner reactions. -^^^ 
It is interesting to note that the CBBN/SBBN ratio of the S'-factor is nearly 0.6 — 1 
in the Gamow peak region in both figures. Also, the ratio was nearly 1 — 2 in the 
previous case shown in Fig. 8. This suggests that the CBBN reaction cross section 
may be approximated, within an error of factor nearly 2, by a simple model in which 
the CBBN S'-factor in (2-21) is replaced by the observed S-factor of the SBBN 
partner reaction but the Coulomb barrier penetration factor is kept the same as that 
in CBBN. 



> 

0) 



LU 



fLiX-) + p 


1 1 1 1 
a + ^He + X" 


SBBN 

-oooooooooooo 


° ° O tl n . 




CBBN " ' 




s-wave 


^V.' It 

,[' 1 J 1 


p-waye 



100 

E [keV] 



200 



Fig. 9. The calculated S'-factor for the CBBN reaction (4-1) (solid line). The s- and p-wave con- 
tributions are shown individually. The observed S'-factor of the SBBN partner reaction (4-3) is 
given by open circles (the adopted values in Ref. 18)). The dotted curve illustrates the Gamow 
peak (in arbitrary units) for Tg = 0.1 {kT = 8.6 keV) with the maximum at Eo = 42 keV, and 
the dot-dashed curve is that for Tg = 0.3 {kT = 26 keV) with Eq = 87 keV. 



This simple model may be used in BBN network calculations (and has already 
been employed in some calculations in the literature) when precise CBBN reaction 
rate is not available. However, the definition of the SBBN partner reaction should 
be applied strictly as discussed in §3.5. 

The S-factor of the CBBN reaction (4-1) in Fig. 9 may be simulated using 
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Fig. 10. The calculated 5- factor for the CBBN reaction (4-2) (solid line). The s- and p-wave 

contributions are shown individually. The observed S'-factor of the SBBN partner reaction (4-4) 
is given by open circles (the adopted values in Ref 18)). The dotted curve is the Gamow peak 
(in arbitrary units) for Tg =0.1 (fcT = 8.6 keV) with the inaxiinurn Eq = 42 keV, and the 
dot-dashed curve is that for T9 = 0.3 (feT = 26 keV) with Eq = 87 keV. 



expression (2-27) with S{E) ~ 2.6 x 10^ keV b. Using (2-28), the CBBN reaction 
rate for Tg < 0.5 is then given as 

Na{(tv) = 2.6 X 10^°r9"^exp(-6.74T9"^) cm^s-^or^ (4-13) 

In Fig. 10, which shows reaction (4-2), the energy dependence of S{E) may be 
approximated by (2-27) with S'(O) = 36 keV b and a = 0.15 b, and therefore the 
reaction rate for Tg < 0.5 is written as 

2 12 

Na{(7v) = 3.5 X 10^ Tg^s exp (-6.74 rg"^) (1 + 0.81 + 0.30 Tg) cm^ s'^ moP^ 

(4-14) 

It will be interesting to see whether, in the BBN network calculation using the 
above reaction rates, the CBBN reactions (4-1) and (4-2) destroy significant amount 
of ^'^Li, since their SBBN partners (4-3) and (4-4) are the reactions that destroy 
^'^Li the most strongly in standard BBN. 

§5. X "-catalyzed charge-exchange reactions 

When the cosmic temperature cools to Tg ~ 0.01, the X~ particle begins to 
form the neutral bound state*-* (pX^) if tx ^ 10^ s; this is also true for (dX^) and 
{tX~), even though their fractions are very much smaller (cf., for example. Fig. 2 in 
Ref. 12)). Owing to the lack of a Coulomb barrier, the bound states are expected to 
interact strongly with other elements that have already been synthesized. 



The binding energy of {pX ) is ~ 0.025 MeV and the rms radius is ~ 50 fm. 
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In this section, we perform a fully quantum three-body calculation of the fol- 
lowing CBBN reactions that are induced by the neutral bound states {pX~), {dX~) 
and (tX-): 

a) Charge-exchange reactions, 

{pX~) + a^ {aX~)3s,3p,3d + P (^'l) 
{dX-)+a^ {aX-)2s,2p + d (5-2) 
(tX-) {aX-)2s,2p+ t (5-3) 

b) Reactions that produce ^'^Li, 

{dX-) + a^ ^Li + X- + 1.4MeV, (5-4) 
{tX-) + a ^ '^Li + X- + 2.4 MeV , (5-5) 

c) Reactions that destroy ^'^Li and ^Be, 

(pX-)+^Li^a+^He+X-+4.0MeV, (5-6) 
{pX-)+'^U^ a+a+X- + 17MeV , (5-7) 
(pX-)+^Be^ ^B+X-+0.11MeV. (5-8) 

In the charge-exchange reactions, the excited states of the {aX") atom with the 
principal quantum numbers n = 3 and 2 are respectively employed in (5-1) and 
(5-2)— (5-3). They have the minimum transferred energy among the energetically 
possible {aX~)ni states. 

Upon performing a DWBA calculation, Jedamzik^^' ^'^^ asserted that the cross 
sections of the CBBN reactions (5-6) and (5-8) that destroy ^Li and '''Be, respectively, 
are so large that (pX^) could induce a second round of BBN, say late-time BBN, 
capable of destroying most of the previously synthesized ^Li and ^Be. However, 
since the strength of the nuclear interactions (on the order of ~ 10 MeV) that cause 
reactions (5-6) and (5-8) is much larger than the incident energy (< 1 keV), the Born 
approximation should not be applied to these low-energy nuclear reactions. 

Since the neutral bound state in the entrance channel behaves like a neutron 
at a large distance with no Coulomb barrier, reactions (5-4) — (5-8) caused by the 
nuclear interaction are expected to have greatly enhanced cross sections even at low 
energies (< 1 keV). Therefore, the most interesting issue in this section is whether or 
not the charge-exchange reactions are so strong that they can intercept the reactions 
(5-4) — (5-8) that produce and destroy Li and Be. 

Before performing a precise quantum three-body calculation, we make a simple 
consideration, following the semiclassical approach in Refs. 14) and 41), to com- 
pare the cross section of the atomic process (5-1) with that of the nuclear reactions 
(5-4) — (5-8). In an analogy to the charge-exchange reactions seen in atomic physics, 
such as the capture of a muon by a hydrogen atom^^^ 

{pe~)is + ^ {piJ'~)nl + e~, (n ~ 14) (5-9) 

*^ The energies of the Is state of {dX~) and (tX~) are —25, —48 and —71 keV, respec- 

tively, whereas those of {aX-)„i are -337 keV(ls),-96 keV(2p),-90 keV(2s),-43 keV(3d),-43 
keV(3p) and -41 keV(3s). 
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we understand that, if the incident a particle in (5-1) enters inside the proton Is 
orbit of {pX~), the sum of the inner charges seen by the proton becomes positive 
(+e), and therefore the proton immediately escapes from the orbit and instead the 
a particle is trapped by as long as the incident energy is not too high. In the 
limit of the semiclassical picture, the cross section of (5-1) is approximately Trfe^, 
where b is the rms radius of the {pX~) atom (6 = 50.4 fm). On the other hand, 
the cross section of the nuclear reaction is roughly 7r6g, where bo is the sum of the 
radius of the incoming nucleus and the nuclear interaction range (bo ~ several fm) . 
Therefore, the atomic cross section is roughly two orders of magnitude larger than 
the nuclear cross section. Actually, this picture is not fully applicable because the 
quantum effect is large for low-lying orbits with n < 3. However, we find that 
the following quantum thrcc-body calculations give a roughly similar ratio of the 
two cross sections, although each shows a strong energy dependence contrary to the 
above simple consideration. 

a) Charge-exchange reactions 

The method of calculating the cross sections of the charge-exchange reactions 
(5-1) — (5-3) is almost the same as that using Eqs. (2-8) — (2-20). For instance, for 
the reaction (5-1), the entrance channel {pX~) + a is considered using the Jacobi 
coordinate system with c = 1 in Fig. 1, whereas the exit channel {aX'~)ni -|- p is 
considered using the system with c = 3. For the closed-channel amplitude ^I'j^^^^'^^ 
in (2-14), we take the same method as that in §2.4. The Coulomb potentials are 
constructed by assuming Gaussian charge distributions of p, d, t and a, the same as 
in the previous sections. Nuclear interactions play a negligible role in these atomic 
processes. 

b) Ractions that produce Li 

The a-transfer reactions (5-4) and (5-5) have the same structure as (2-2) and 
(2-4), respectively. Therefore, the method for calculating the cross sections of the 
former reaction is the same as that in §2 except that the entrance channel here is 
described using the Jacobi coordinate system with c = 3 in Fig. 1. 

c) Reactions that destroy Li and Be 

The two reactions (5-6) and (5-7) have the same structure as (4-3) and (4-4), 
respectively. Thus, the method for calculating them is the same. The structure of 
reaction (5-8) is similar to that of (3-2), although and X~ are not in a bound 
state here but in a photonless scattering state. In (3-2), the transition from the 
entrance channel to the exit channel via the electromagnetic interaction was treated 
perturbatively, but reaction (5-8) is solved as a ^Be + p + X~ three-body problem. 

5.1. Results 

The calculated cross sections cr(E) of reactions (5-1) — (5-8) are listed in Table II 
for energies E =0.01, 0.1, 1 and 10 keV (Tg - 0.0001,0.001,0.01 and 0.1, respec- 
tively). Since the Coulomb barrier penetration factor exp{—2TTrj(E)) in (2-21) is 
unity here, the astrophysical S'-factor is simply given by S{E) = Ea{E). 

The extreme enhancement of the cross sections of reactions (5-4) — (5-8) due to 
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the neutral bound states can be seen, for example, by the fact that the Li production 
cross section of (5-4) at = 10 keV is 6.4x 10~^ b, whereas that of (2-2) with {aX~) 
in the entrance channel is 3.9 x 10~^ b.^^ As can be seen in the usual neutron-induced 
low-energy reactions, the cross sections roughly follows the 1/v law at lower energies 
( < 1 keV) in Table II. 



Table 11. Calculated cross sections and reaction rates of the late-time BBN reactions induced by 
the neutral bound states. The rates (in units of cm^ mol~^) are applicable for Tg ^ 0.05. 







cross 


section (b) 
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Reaction 


0.01 kc\' 
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a) charge-exchange reaction 












{pX-)+a^ {aX-):u+p 


8.4 X 10^ 


2.2 X 10^ 


7.8 X 10^ 


7.5 X 10^ 


1.0 X 10^° 


idX-) + a^ {aX-}2e + d 


3.1 X 10^ 


9.1 X 10^ 


2.0 X 10^ 


2.1 X 10^ 


3.5 X 10^ 


(tX-) + {aX-)2e+ t 


8.3 X 10^ 


1.9 X 10^ 


3.2 X 10^ 


2.5 X 10^ 


7.6 X 10^ 


b) a-transfer reaction 












{dX-) + a^ ^Li-I-X- 


9.6 X 10° 


3.0 X 10° 


6.9 X 10"^ 


6.4 X 10"^ 


1.1 X lO'' 


(tX-) + '^U + X- 


3.5 X 10"^ 


1.1 X 10"^ 


2.7 X 10"^ 


3.0 X 10"^ 


4.3 X 10^ 


c) Li-Be destruction 












(pX-)-h®Li^ a-l-^He-l-X- 


1.8 X 10^ 


5.5 X 10^ 


1.1 X 10^ 


2.8 X 10° 


1.6 X 10® 


{pX~)+'^U^ a+a+X~ 


3.8 X 10° 


1.7 X 10° 


7.7 X 10"^ 


1.6 X 10"^ 


5.5 X 10® 


(pX )+^Be^ ^B+X- 


4.8 X 10"^ 


3.2 X 10° 


5.3 X 10"^ 


5.5 X 10"^ 


5.2 X 10® 



The most important result in Table II is that the cross sections of the charge- 
exchange reactions are much larger than those of the nuclear reactions (5-4) — (5-8). 
Therefore, for the a-transfer reactions (5-4) and (5-5), the bound states {dX~) and 
{tX~) are mostly changed to {aX~) before producing ^'^Li. For the Li-Be destruc- 
tion reactions (5-6)— (5-8), we further note that the probability that {pX~) comes 
in contact with ^'^Li and ^Be is more than several orders of magnitude smaller than 
that of it coming in contact with an a particle owing to the large difference in the 
element abundances, and that, even if (pX)^ did collide with an element A (=^'^Li 
or ^Be ), another atomic charge-exchange reaction, 

{pX-) + A^{AX-)ni+p, = ^•^Li,^Be) (5-10) 

would immediately occur instead of reactions (5-6) — (5-8) since the cross section of 
(5-10) would be as large as that in (5-1), although it was not calculated here. 

To conclude, the nuclear reactions (5-4)— (5-8) that produce and destroy Li and 
Be would negligibly change the element abundance in the late-time BBN owing to 
the strong interception by the charge-exchange reactions (5-1) — (5-3) and (5-10). 

We derive the reaction rates of (5-1) — (5-8) for use in the network BBN calcu- 
lation under the approximation that the cross sections in Table II may be roughly 
simulated in a simple form (the 1/v law) 



(5-11) 
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An evident reduction of the cross section in Table II from (5-11) can be seen at -B ~ 10 
keV due to the structure of the neutral bound state and the decreased probability 
of the a particle being trapped by X~ , but this reduction is thought to have little 
effect on the clement abundance*-* . The use of this approximation in the calculation 
of the reaction rate (2-23) makes the result independent of the temperature Tg.**^ 
The rate for Tg < 0.05 is described as 

NA{crv) = 2.6 X lOV^m//^ (t{E^) cm-'^s^^mol'S (5-12) 

where the constant C is determined by matching (5-11) to cr{E) in Table II at 
E = E^. Here, E^ and cr(£'m) are expressed in units of keV and b, respectively, 
and jU is the reduced mass of the entrance channel in units of amu. We then simply 
averaged the three rates determined at = 0.01,0.1 and 1 (keV). 

The calculated reaction rates are listed in the last column of Table II. We consider 
that it is not necessary to simulate a{E) in a more sophisticated manner than (5-11), 
because it is clear from Table II and the above discussion that the late-time BBN 
does not affect the element abundances meaningfully***^. A similar conclusion on 
the effect of the late-time BBN was obtained in Ref. 14). 

§6. X "-catalyzed production of ^Be 

The absence of stable isotopes with mass number 8 is a bottleneck for the produc- 
tion of ^Be and heavier nuclei in standard BBN. In recent papers, Pospelov and 
coworkers pointed out the possibility of the enormous enhancement of the reaction 
rate for the production of ^Be by considering an X~-catalyzed resonant neutron- 
capture process 

(^BeX") + n ^ CBe,+X~) ^Bcg- + X' . (6-1) 

2 2 

They claimed that the primordial abundance of ^Be imposes strong restrictions on 
the lifetime and abundance of X~ . Here, ^Bca- is the ground state and ^Bei + 

2 2 

denotes the first excited state with spin ^ at Ex = 1.735 it 0.003 MeV with the 
neutron decay width r„ = 0.225 ± 0.012 MeV^^) (or E^ = 1.684 ± 0.007 MeV with 
r„ = 0.225 ± 0.012 MeV^^)). 

In this section, we make a critical comment on this work.-'-^)'-'-^) Their attention 
to reaction (6-1) is interesting, but there is a serious problem in their calculation. 

*^ Also, the cross section at E = 0.01 keV in the bottom Une of Table II is much smaller than 
that given by (511), but this is because the p-wave contribution is dominant in (5-8) and momentum 
matching becomes difficult at low energies when the proton is transferred to a definite bound state 
in the exit channel. 

**' Pospelov et o?.^*' derived the rate as oc l/v/Tg, but this is due to the fact that they did not 

calculate the energy dependence of the cross section and assumed it to be constant with respect to 
the energy. Their value at Tg = 0.01 is of the same order of magnitude as ours for {pX~) + a —^ 

(aX~)exc. +p- 

***•* According to a BBN network calculation by Kusakabe^'^' including the reaction rates in Table 
II, the abundances of the bound states {pX~), {dX~) and {tX~) at Tg ^ 0.05 are greatly reduced 
by a factor of ~ 10® from those obtained without including the reaction rates. 
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They assumed that the rms charge radius of both and ^Be was 2.50 fm. Since 
the observed charge radius of ^Bes- is 2.519 ± 0.012 fm,^) the assumed radius is 

2 

acceptable for ^Bes-. However, there is no reason to take a radius of 2.50 fm for 

2 

^Be, which is a resonance state at 0.092 MeV above the a + a threshold. In principle, 
the radius of such a resonance state cannot be experimentally measured, and even 
theoretically, it is difficult to define the radius of the oscillating resonance wave 
function. Since the width of this resonance is, however, very small (5.57 eV), the wave 
function is heavily attenuated by the Coulomb barrier, followed by an asymptotically 
oscillating amplitude that is roughly three orders of magnitude smaller than that in 
the nuclear interaction region. Therefore, it is not meaningless to derive the rms 
radius by using the resonance wave function with the asymptotic part omitted, or 
by using the wave function obtained through the diagonalization of the Hamiltonian 
with appropriate L^-integrable basis functions; both methods result in almost the 
same radius as long as they give the same resonance energy. 



1.168 MeV Be + n + X 

^ 

: 0.092 

1 .076 a + a + n + X 



1.574 



0.0 MeV 



("Be X") + n 



-0.233 (-0.284) 
(^Be,/2,X-) I 



1.735 
(1 .684) 

-1.915 



-0.445 



t "Be3,2_ + X- 
1.470 



("Be3,2_ X") 

Fig. 11. Isobar diagram for the *Be+n+X~ system. This is to be compared with Fig. 1 of Ref. 13); 
here the same model as that in Ref. 13) is employed but the rms charge radius of *Be is changed 
to an appropriate value from the too small radius used in Ref. 13) (see text). It is demonstrated 
that the (^Bei + X~) state becomes below the (*BeX~) -|- n threshold contrary to Ref. 13). 



A microscopic a+a model calculation using the resonating group method (RGM) 
by Aral et al.^^^ gives 3.16 fm for the rms charge radius of ^Be but the resonance is 
0.022 MeV above the a + a threshold. A similar calculation with the orthogonality- 
condition model (OCM)^^) by the present authors and their coworkers^^^ gives a 
radius of 3.39 fm when the a + a resonance is at the observed position of 0.092 MeV. 
Since the latter model gives 3.20 fm when the resonance energy is tuned to the same 
value (0.022 MeV) as that in the former RGM calculation, we here employ an rms 
charge radius of ^Be of 3.39 fm. 

We now repeat the same analysis as that in Refs. 13) and 14), using the same 
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model but with charge radh of 3.39 and 2.52 fm for ^Be and ^Bc, respectively. Here, 
mx — oo is taken similarly to Ref. 13). As illustrated in Fig. 11, the binding 
energies of (^BeX~) and (^Bes-X") are respectively 1.168 and 1.470 MeV with 

2 

respect to the ^Be + X" threshold and ^Bes- + X~ threshold that is located 1.574 

2 

McV below the a + a + n + X~ threshold. If we assume, similarly to in Ref. 13), 
that '^Bci+ has the same charge radius as ^Bea-, then the (^Bei + X~) state is 1.735 

2 "22" 

MeV^3) (oj. i_gg4 MeV^"^)) above {^Be^-X-). As a result, we find that (^Bei+X") 

2 2 
is not a resonance state but a bound state located 0.233 MeV (0.284 MeV) below 

the (^BeX~) + n threshold with no neutron decay width (although the state has 

a small width allowing decay to the ^Bes- + X~ channel). Those level energies 

2 

depend on the shape of the charge distributions of ^Be and ^Be, but the effect on 

the energy difference between the (^Bei+X~) state and the (^BeX~) +n threshold 

2 

is very small; for example, the difference changes by only 0.005 MeV if the charge 
distribution is artificially changed from Gaussian to square in both the nuclei. We 
thus conclude that the resonant neutron-capture process (6-1) is not realized as long 
as the model in Ref. 13) is employed. 

Another problem in the analysis in Refs. 13) and 14) is that the charge radius of 

the excited resonance state ^Bei+ is assumed to be the same as that of the ground 

2 

state of ^Be3-. An a + a + n three-body RGM calculation by Aral et al^^^ gave 2.88 
2 

fm for the rms charge radius of ^Bei + . The calculated energy is ~ 0.4 MeV above 

2 

the a + a+n threshold while the observed value is 0.161 MeV^^^ (0.110 MeV^^)). The 
calculated neutron-decay width of ~ 0.2 McV almost reproduces the experimental 
value. According to Aral et a/.,^^-* an artificial change of the resonance position, even 
down to slightly below the a + a + n threshold (by tuning a part of the interaction), 
has little effect on the charge radius of ^Be 1 + ; this is because the change in the wave 

2_ 

function is dominantly seen in the resonating valence neutron and little change is 
seen in the a clusters. If we employ the charge radius of 2.88 fm for ''Bci + , the 

energy -0.233 McV (-0.284 McV) of (^Bci+X") shown in Fig. 11 must be replaced 

2 

by —0.116 MeV (—0.167 MeV). Therefore, the consideration of the differences in the 

rms charge radii among ^Be, ^Bes- and ^Bei+ still causes the (^Bei+X~) state to 

22 2 
lie below the a + a + n threshold. 

In the above discussions on (6-1), it is still necessary to consider the change in 

structure of the loosely coupled systems ^Be and ^Bei+ induced by the addition of 

2 

X~ to them. The dynamics can only be studied by employing an a + a + n + X~ 
four-body model.*) Use of this model will also make it possible to exactly calculate 
the (^BeX~) formation processes a + X~ — {aX~) + j, {aX~) + a — (^BeX~)+7 



*' This type of structure change by the injection of an impurity particle into loosely coupled 
nuclear states has well been studied in light hypernuclei by the present authors and their cowork- 
ers'*®^"*®^ using three- and four-body models. 
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and the resonant (nonresonant) ^Bes- formation proeess (6-1). This study based on 

2 

the four-body model is beyond the scope of the present paper, but such a study is 
in progress and will be presented in the near future. 

§7. Summary 

(1) Using a fully quantum three-body method developed by the authors/^) we 
have calculated the cross sections of various types of big-bang nucleosynthesis (BBN) 
reactions that are catalyzed by a hypothetical long-lived negatively charged, massive 
leptonic particle (called X~) such as the supersymmetric (SUSY) particle stau. We 
also provided their reaction rates for use in the BBN network calculation. The rates 
arc summarized in Table III and in the last column of Table II (the late-time BBN 
reactions for Tg ^ 0.05). 

(2) In the catalyzed BBN (CBBN) reactions, the strength of the nuclear interac- 
tion (on the order of 10 MeV) , which causes the transition between the entrance and 
exit channels, is very much larger than the incident energy (< 200 keV), and there- 
fore the coupling between the two channels is sufficiently strong to induce multistep 
transitions between the channels. Therefore, any calculation method that does not 
take into account the above property of the interaction is not suitable for application 
to the reactions. Also, for the charge-exchange reactions (5-1)— (5-3), a fully quan- 
tum, nonperturbative treatment is highly desirable since the a-particle is transferred 
to low-lying states with the principal quantum number n = 2 — 3. We have shown 
that our three-body calculation method^^^ satisfies the above requirements and is 
useful for the study of all the reactions. 



Table III. Summary of the calculated reaction rates of CBBN reactions obtained by the three-body 
calculation, a) - c) are for Tg <, 0.2 and d) - g) are for Tg 0.5. 



Reaction 



Reaction rate (cm'^ s ^ mol ^) 
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(3) In the typical photonless CBBN reactions a)— c) in Table III, the calculated 
S'-factors have similar orders of magnitude to those of the typical nonresonant pho- 
tonless reactions in standard BBN (SBBN). As was first pointed out by Pospelov^^ 
and then confirmed by Hamaguchi et al,^^ the S'-factor of CBBN reaction a) that 
produces ^Li is many orders of magnitude larger than that of the radiative capture 
a + d — ^^Li + 7 in SBBN. This imposes strong restrictions on the lifetime and abun- 
dance oi X~ . However, this large enhancement is simply because this SBBN reaction 
is heavily El-hindercd. On the other hand, the S'-factors of CBBN reactions b) and 
c) are found to be only ~ 30 times larger than those of the strong El radiative 
capture SBBN reactions a + t(^He) — >-^Li(^Be) -|- 7. Therefore, reactions b) and c) 
do not seem to change the abundances of ^Li-^Be meaningfully. 

(4) The resonant radiative capture CBBN reaction, g) in Table III (cf. Fig. 5), 
was proposed by Bird et al.^^ as a possible solution to the overproduction of ^Li-^Be. 
We have examined their model and result using a ^Be + p + X~ three-body model, 
which makes it possible to calculate all the steps in the reaction. We have shown 
that both the very loosely coupled ^Be + p structure of ^B and the strongly spin- 
dependent ^Bes- — p interaction contribute largely to the energy of the resonance 

2 

{^BX^Y^p- However, the effects work oppositely, almost cancelling each other, and 

therefore the reaction rate obtained by the simple model^^ (assuming a Gaussian 

charge distribution of ^B) is by chance close to ours. The reaction rate is sensitive 

to nix', the rate with mx = 100 GeV is 60% of that with nix ^ 00 at Tg = 0.3. 

Regarding the recombination process that forms (^Be-^-X^) starting from ^Bes- 

2 2 

and a proton, we pointed out in §3.6 that the calculation^^ by Bird et al. of the energy 

of the 'resonance' state (^Bei-X~)2s is erroneous owing to the lack of consideration 

2 

of the nuclear structure of ^Bca- and ^Bci-; the a + ^He + X~ three-body cal- 

2 2 

culation clarifcd that the (^Bei-X~)2s state cannot be a resonance but must be a 

2 

bound state below the ^Bcs- + X~ threshold. In Ref. 6) two types of results were 

2 

reported on the abundance of '^Li-^Be calculated with and without the (^Bei-X^)2s 

2 

resonance, but the case with the resonance is not acceptable. 

(5) We have performed fully quantum mechanical study on the late-time BBN 
reactions (5-1)— (5-8), which take place at Tg < 0.01 between the neutral bound 
states, {dX~) and (tX~)^ and light nuclei. The atomic charge-exchange 
reactions (5-1) — (5-3) are so dominant that they immediately intercept the nuclear 
reactions (5-4) — (5-8). Therefore, we conclude that change of element abundances 
due to the late-time BBN reactions is negligible. 

(6) We have examined the 'resonant' neutron-capture CBBN process (6-1) for 

producing ^Be proposed by Pospelov and coworkers. ^^^'^^^ They claimed that the 
resulting increased output of '^Bc strongly constrains the abundance and lifetime 
of X^ . However, we found that, as long as the same model as that in Ref. 13) is 
employed but with an appropriate charge radius of ^Be (a too small value of the 
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radius was used in Ref. 13)), the intermediate state {^Bci + X ) appears below the 

2 

(*BeX~) +ra threshold, thus giving no resonant mechanism. The model is, however, 
too simple to account for the change in the structure of the loosely coupled (resonant) 
systems ^Be and ^Bei+ due to the injection of the X" particle into them. Therefore, 

2 

an extended four-body calculation based on an a + a + n + X~ model is underway 
to examine (6-1) much more seriously. 

(7) Finally, we expect that the application of the presently obtained reaction 
rates to the BBN network calculation will help solve the ^Li-^Li problem and simul- 
taneously impose restrictions on the primordial abundance and lifetime of the X" 
particle. 
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